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Objectives  of  the  Project 

Plasma  chemistry  mechanisms  for  accelerating  activation  of  fuels  by  projected  plasma  jets 
are  investigated  by  the  team  of  Moscow  State  University  (MSU).  Main  attention  is  paid  to  excited 
state  transitions  and  related  dissociation  processes  inside  the  excited  area  created  by  a  plasma  jet 
and  in  its  vicinity.  This  approach  is  different  from  the  conventional  paradigm  that  uses  plasma- 
dynamic  discharge,  stimulating  fuel  ignition  with  short  powerful  ultra-violet  radiation.  It  is 
envisioned  that  this  new  approach  will  produce  substantially  more  effective  ignition  than  the 
conventional  method  at  much  lower  temperatures,  and  without  decomposition  of  the  fuel  by 
radiation,  which  take  place  in  the  known  method.. 

For  verification  of  the  proposed  hypothesis  we  will  undertake  plasma  jet  parameter 
experimental  investigations  in  the  motionless  is  supposed,  theoretical  and  computation 
investigations  of  air-propane  mixture  at  pressures  300-760  with  a  help  of  nitrogen  plasma  jets  are  to 
be  executed.  The  effect  of  plasma  activation  on  the  ignition  delay  time  and  combustion  velocity  will 
be  theoretically  investigated. 

The  objective  of  this  project  is  to  establish  the  feasibility  of  plasma-jets  for  plasma  chemical 
stimulation  of  combustion  and  aerodynamic  augmentation.  Particular  attention  will  be  focused  to 
the  understanding  of  the  plasma  chemical  mechanisms  of  the  activation  of  fuels,  coupling  with  the 
fluid  dynamic  processes  and  their  implications  for  aerodynamic  and  propulsive  augmentation. 
Specific  problems  to  be  solved  are: 

•  Assessments  of  new  scheme  for  using  atoms  and  active  radicals  for  fuel  activation 

•  New  experimental  data  necessary  for  oblique  plasma  jet  in  flow  application  for  reinforcement  of 
new  scheme  from  Item  1 . 

•  Results  of  numerical  modeling  interaction  of  oblique  plasma  jet  with  streams  with  and  without 
fuel  mixing  and  burning. 

•  Recommendations  definitions  of  high  speed  fuel  flows  ignition  and  combustion  stabilizing 
conditions  optimization  on  a  basis  of  numerical  modeling  and  experimental  data  on  plasma 
kinetic  processes  from  Item  1  (in  combination  with  plasma  jets  oblique-to-the-free  stream 
injection  and  flow  gasdynamics  ). 
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Introduction 

According  to  the  concept  of  om'  proposal  plasma  can  be  applied  for  improvement  of 
combustion  processes  in  air  mixbires  with  organic  components  (propane),  in  so  doing  it  excites  a 
gas  flowing  in  a  combustor.  Plasma  generators  are  applied  for  plasma  creation.  Ionization  and 
dissociation  of  molecules  takes  place  in  plasmas.  Active  paificles-active  radicals  and  atoms- 
products  of  dissociation  of  molecules  are  present  in  plasmas. 

Mixing  of  plasma  particles  with  fiiel  gas  molecules  takes  place  at  interaction  of  plasma  jet 
with  mixture  flow  takes  place  in  a  combustor.  These  particles  posses  high  reaction-ability  and  they 
play  important  role  in  chemical  reactions  quickening  oxidation  process  at  presence  of  oxygen 
molecules  and  atoms.  At  presence  of  hydrocarbon  components  in  a  mixture  such  processes  lead  to 


Fig.  1.1  Principle  scheme  of  plasma  fuel  activation  with  a  help  of  plasma  generator. 


(«insulator»),  6  -  plasma  generator,  7  -  combustor  (section  JSfo2),  8  -  measmement  section 

In  Figur  e  1.1  one  can  see  a  principle  scheme  of  plasma  generator  application  for  combustion 
activation  in  a  gas  flow.  In  Figure  1.2  one  can  see  a  supposed  experimental  scheme  in  which 
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proposed  concept  can  be  realized.  In  this  figure  one  can  see  a  test  channel  for  investigations  of 
plasma  jet  in  a  supersonic  flow,  where:  1.  is  a  duct  for  creation  of  required  flow  characteristics;  2.  Is 
storage  and  supply  air  high  pressure  system  at  pt,  max  ^  12;  3  -storage  and  supply  gaseous  fuel 
(propane)  system;  4.  plasma  generator  and  system  of  electrical  energy  input;5.  measurement  system 
of  gasdynamic  and  electrical  parameters;  6.  system  of  visual  observation  and  video  taping  of  flow 
visualization  in  test  channel. 

Development  and  manufacturing  of  highly  effective  technological  plasma  generators  for  plasma 
chemical  activation  of  supersonic  and  subsonic  fuel  flows  and  obtaining  high  enthalpy  plasma  flows 
with  their  help  leads  to  necessity  of  consideration  plasma  generator  design  with  vortex  flow 
stabilization  and  divergent  anode  duct,  see  Fig.3.  Most  important  features  of  such  geometry  are  high 
consumption  characteristics,  efficiency  of  working  gas  heating  and  smallness  of  heat  losses  into 
water  cooled  anode  surface  [1.1]. 

Plasma  generators  of  this  type  differ  from  others  by  high  enthalpy  plasma  jet  generation 
capability  with  high  consumption  of  a  working  gas  at  insignificant  internal  jet  length  (10-40  mm) 
and  high  coefficient  of  efficiency.  Since  maximum  temperature  value  in  such  generators  with 
typical  diameter  of  the  inletting  anode  hole  3-5  mm  can  exceed  20000  K,  and  average  mass 
temperature  at  the  outlet  reaches  4000-6000  K  at  a  gas  consumption  0,5-  5  g/s,  then  these  plasma 
generators  have  wide  technological  applications  [1.1-4]. 

Reliability,  comparable  simplicity  of  design  of  these  plasma  generators  and  considerable 
amount  of  experimental  data  on  plasma  characteristics  starting  from  near  cathode  region  of  electric 
discharge  and  finishing  with  region  of  plasma  relaxation  at  a  distance  up  to  50  from  the  anode 
down  a  flow.  This  became  convincing  bases  for  choice  of  mainly  this  scheme  of  plasma  generator’s 
design.  High  efficiency  of  electric  field  energy  in  the  arc  to  thermal  energy  of  a  flow  is  insured  by 
vortex  jet  stabilization  at  tangential  delivery  of  plasmaforming  gas  (nitrogen,  argon).  Convective 
and  conductive  loses  to  water  cooling  wall  of  diverging  duct  are  considerably  decreased  at  such 
means  of  the  discharge  organization.  So  main  loses  are  connected  with  radiation  energy  losses  and 
their  role  increases  with  rising  of  plasma  temperature. 

Most  convenient  plasma  forming  gases  from  the  point  of  view  of  repeated  use  of  plasma 
generators  are  nitrogen  and  argon.  They  insignificantly  destroy  electrodes  in  comparison  with  air. 

Calculations  of  thermal  plasma  composition  of  air  and  nitrogen  of  ours  our  and  [1.1-7],  see 
Fig.  1.4-5  allow  to  conclude  that  main  prospective  temperature  range  at  plasma  generator’s  outlet  is 
T=  4000-10000  K  from  the  point  of  view  active  radical  generation.  Calculations  show  that  atom 
concentration  is  of  the  order  of  magnitude  of  molecules  or  higher  in  this  temperature  range,  so  one 
have  expect  here  the  most  efficiency  of  realization  of  the  proposed  concept. 
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Availal  ^  -  Geometiy  of  electrode  units  of  plasma  generator:  1  -  cathode,  2  -  water 
cooled  anode  with  diverging  duct,  3  -  plasma  forming  gas  (Ni?  Ai). 
in  range  r-iu  -  u  kdt. 


.  ,  Fig.  1.5  Component  composition  of 

Fig.  1,4  Component  composition  of  equilibrium  ail' plasma  at  constant  pressure 

equilibnum  nitrogen  plasma  at  constant  p-lO^P 

pressuie  P=10^  Pa 


Chapter  1.1.  Main  arc  plasma  equations  applied  in  calculations 
At  considerably  high  pressuie  (about  of  the  atmospheric  one  and  higher)  electiic-arc  is  close  to  the 
thennal  equilibrium  one  [1.4, 1.7],  The  evaluation  of  the  deviation  of  electron  temperatuie  Tg  from 
the  temperature  of  heavy  pailicles  Tg  can  be  made  by  known  relation  [1.7-8] 
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)  (1-1) 

Te  ~  32  3^^ 

U  V 

where  wjg,  Wg,  are  masses  of  a  gas  atom  and  electron,  e  is  electron  charge,  Xg  is  electron  mean  free 
path,  E  -  electric  field,  k  -  the  Boltzmann  constant.  It  can  be  transformed  to  more  obvious  form,  if 
we  consider  that  the  Joule  heat  from  the  arc  is  transported  mainly  by  the  thermal  conductivity, 
T-crE^r^  IX  , 

and  take  into  the  account  the  equation  for  the  current 
/  ^cETir^ , 

where  cris  plasma  conductivity,  X  -  thermal  conductivity  of  a  gas,  r  is  characteristic  radial  size  of  an 
arc 

(  V 

4^-^  .  (1.2) 

T’.  ~32m,  3^^ 

(  2  j 

For  typical  values  of  nitrogen  at  the  temperature  Tg  -2000-  4000  K,  X  -0.1  W/(m  K),  Xg  -10'^  m 

T 

already  at  /  -  10  A  we  have  1 — ^-10"^,  or  plasma  state  within  the  arc  is  close  to  thermal 
equilibrium. 

Because  of  this  result  one-temperature  model  of  plasma  in  arc-discharges  is  usually 
developed  and  we  also  will  use  it  in  our  approach.  For  this  situation  at  the  first  step  of  analysis  we 
will  represent  the  equation  system. 

Gasdynamic  part  of  equation  system  for  electric-arc  discharge  plasma  includes  the  continuity 
equation 

^  +divip-V)=0;  (1.3) 

the  momentum  equation  is 

PV  -  -  ^  ^  ^  ^ 

P  —  +  p{V-grad)V=ip-pJg+EdivD+jxB- 
d  t  ;  (1.4) 

-grad{P+  (213)7]  divV)+2div  (t]  S) 
the  energy  equation  is 
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(1.5) 

(1.6) 

(1.7) 


divB  =  0,  (1.8) 

and  Ohm’s  law 

E  +  V  xB  =  ^  +  —  y]  ><B-gradp\  (1.9) 

cr  en^ 

System  (1.3)-(  1.9)  is  supplemented  with  the  equations 
p^p{T,P)  ;  a  ^cj{T,P)  ;  1  ^X{T,P)  ;  rj  (T,P)  ;  (1.10) 

=c^  (T,P)  ;  h  =h  (T,P)  ;  'F  ='F  (T,P)  ;  (1.11) 

B  =  p^H;  D  =  £^E.  (1.12) 


The  following  notations  in  (1.3)-(  1.12)  were  used:  g  is  a  vector  of  the  free  fall  acceleration;  S  is 
the  tensor  of  deformation  velocities  with  the  components 

^ik  =(^K-  /5x,.)/2, 


where  ,  V^.  are  the  components  of  the  vector  V ;  x. ,  are  coordinates,  k,  i  =  l,  2,  3 ;  p^  is 
density  of  the  ambient  medium,  P^  is  partial  density  of  electrons,  'F  is  radiation  energy  flux. 

Let  us  simplify  the  system  of  equations  (1.3) -(1.12).  Estimates  show  that  in  equations  we 
can  neglect  by  the  Archimedes  and  the  Coulomb  forces.  Indeed  the  Archimedes  number 
Ar  =  \p-p^\gLI{pV^) 

{L  is  characteristic  length).  For  typical  conditions  of  arcs,  p  ~5-10”^  kg/m^,  -1.5  kg/m^,  V ~ 
10^  m/s,  L  ~  10'^  m  we  have^r-  10'^. 

In  the  Ohm’s  law  it  is  usually  enough  to  account  the  first  terms  in  the  left  and  right  hand 
sides  of  the  equation.  Indeed,  the  estimates  of  the  typical  conditions  at  P=\  atm,  I  =  200  A,  (cr-  lO"^ 
sm/m,  F~  10  m/s  -  10'  m  ,  ng  -  10  m'  ,  Z  -  10'  m)  show  that  the  ratios  of  density  of  the 
inductive  electric  current,  density  of  Hall’s  current  and  density  of  the  current  caused  by  the  gradient 
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of  the  electron’s  pressure  to  the  density  of  the  current  in  the  arc  have  the  following  orders  of 
magnitude 
gFB 


H,aVL~10-^  ; 


J 


aJB  Upg/ 
jen^  Inenfi 

In  such  cases  the  represented  terms  can  be  omitted  because  of  their  small  value  and  the 
Ohm’s  law  will  be  used  in  its  simplest  form 

]  =  aE.  (1.13) 

Using  (1.7)-(  1.13)  we  can  determine  the  ratio  of  the  electrostatic  force  to  the  electromagnetic  one 


E 

div  D  f 

~  *0 

j 

B 

where  2S  is  the  thickness  of  the  conducting  region  of  the  arc,  ,S~  10'  m. 

Further  simplification  of  the  system  (1.3)-(  1.12)  are  connected  with  the  symmetry  condition 
d !  d(p  =  0 

and  the  condition  of  absence  of  external  magnetic  field. 

Boundary  conditions  for  the  considered  system  of  equations  are  the  following. 

Existence  of  the  attachment  of  the  flow  on  the  walls 

V  =  0,  (1.14) 

mall  thermal  flow  to  walls 


on 


(1.15) 


where  n  is  the  normal  to  the  wall’s  surface. 

Tangential  component  of  electric  field  satisfies  the  following  conditions  on  the  metallic 

walls 

E^=0.  (1.16) 

On  the  dielectric  surfaace  and  between  electrodes 
/ 


=- 


(1.17) 


Where  I  the  current, 
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on  the  outlet. 

On  the  symmetry  axis  we  have 
dP  dT 

V^=0,H^=0,—=0,  — =0. 

or  on 


(1.18) 


(1.19) 


This  equation  system  and  eorresponding  boundary  conditions  are  used  in  our  conditions. 

In  calculations  for  the  Project  we  use  reliable  thermal  physical  data  and  their  interpolations.  Data 
for  density  p(  p,  T),  specific  heat  =c^  (P,P)  and  enthalpy  h{  p,  T)  for  nitrogen  and  air  are 

represented  in  [1.3-5,  1.11-17  ].  Data  about  viscosity  ri(  p,  T),  electric  conductivity 
(7  =(7{T,P)  and  thermal  conductivity  A( p,  T)  are  presented  in  [1.3-5].  Data  on  radiation  'Po  are 
presented  in  [1.3]  their  detailed  analysis  one  can  find  in  [1.3,  1.4,1.17]. 

Analysis  of  experimental  and  theoretical  data  for  nitrogen  showed  that  at  pipn  =  0.3  -  3, 
where  p„  is  normal  atmospheric  pressure  one  at  calculations  can  use  the  following  approximations 
p(p>  T)=p(p„,  T)  pIpn, 

Po(p,  T)  =  Po(Pn,  T)  ■  [p(p,  T)/p(p„,  T)f, 
h{p,T)  =  h(p„,  T),  pi p,  T)  =  p( pn,  T), 

Ap,T)  =  X(p„,  T), 

The  mass  flow  rate  is  given  at  the  inflow  to  the  plasma  generator,  pv^ ,  ,  the  pressure  is 

extrapolated  inside  the  plasma  generator  to  the  inflow  boundary,  ^  =  density  and  the 

dz 

temperature  are  determined  from  the  condition  of  the  isentropic  gas  transition  from  the  normal  state 

to  the  state  with  the  determined  pressure,  «  v  v  (S-  is  the  entropy,  y  -  is  the  adiabatic  exponent), 

/P 

the  velocity  of  the  inflow  is  determined  by  the  given  mass  flow  rate. 

Necessary  data  for  calculations  of  interaction  of  plasma  jets  with  air  and  nitrogen  are 
presented  in  the  following  references. 

Sources  for  thermal  physics  data 


Gas 

P 

N2 

[1.11] 

Air 

[1.13-16] 

ii\  \\f  h 

[1.12]  [1.13]  [1.11] 
[1.4]  [1.16]  [1.14-16] 


X 

Cp 

a 

[1.12] 

[1.11] 

[1.12] 

[1.4] 

[1.14-16] 

[1.4] 
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Chapter  1.2.  Free-Lagrange  Method  for  Computing  of  Two  -  Dimensional  Magneto-Gas- 

Dynamic  Flows 

In  present  investigations  we  use  the  Implicit  Free-Lagrange  method  for  carrying  out  of  our 
Computations.  Let  us  briefly  stop  on  it. 

The  IFLA  (Implicit  Free-Lagrange)  method  was  developed  for  computing  two-dimensional 
non-stationary  MGD  flows  (see,  e.g.  [1.9-10]).  This  is  a  method  based  on  an  implicit  completely 
conservative  operator-difference  MGD  on  an  irregular  triangular  Lagrangian  grid  of  variable 
structure.  The  grid  structure  must  be  changed  because  a  Lagrangian  grid  associated  with  a  moving 
medium  may  be  distorted  in  the  process  of  computation.  The  appearance  of  triangular  cells  that 
strongly  deviate  from  a  regular  triangle  disrupts  the  approximation  of  differential  operators  by  their 
grid  analogs  and  ultimately  leads  to  breakdown  of  the  grid.  Automatic  grid  adjustment  in  computing 
time,  which  involves  arbitrary  changes  of  grid  structure  preserving  the  triangularity  property,  makes 
it  possible  to  maintain  a  nearly  perfect  grid.  Grid-structure  adjustment  in  the  IFLA  method  is 
possible  also  because  the  difference  operators  are  written  on  an  arbitrary  irregular  triangular  grid. 

The  non-stationary  IFLA  method  is  a  sequence  of  time  steps.  The  n-th  step  input  is  an 
irregular  triangular  grid  with  the  values  of  magneto-gas-dynamic  functions  at  the  instant  t„ 
defined  in  its  cells  and  nodes.  Each  time  step  is  executed  in  two  stages.  The  first  stage  analyzes  the 
quality  of  the  triangular  grid,  and  if  necessary  makes  a  local  adjustment  of  grid  structure  on  the 
basis  of  the  metric  principle.  The  function  values  are  recalculated  for  the  new  grid  )’,  so  that  the 
recalculation  for  the  centers  of  the  new  cells  preserves  the  conservation  laws  and  the  gradient  values 
of  these  functions  at  the  nodes.  The  latter  condition  is  highly  significant  for  computation  of 
complex  magneto-gas-dynamic  flows.  The  second  stage  of  the  time  step  evaluates  the  magneto-gas- 
dynamic  functions  on  the  grid  at  the  instant  4+7  using  an  implicit  operator-difference  scheme. 
The  implicit  Completely  Conservative  Operator-Difference  Scheme  (CCDS)  of  two-dimensional 
MGD  used  in  our  method  is  constructed  and  analyzed  by  operator  approach  [1.18]  of  the  theory  of 
difference  schemes,  which  has  been  developed  and  consistently  implemented  in  all  stages  of  our 
work. 

IFLA  method  was  thoroughly  tested  in  solution  of  many  problems  and  proved  to  be 
successful  at  solution  of  important  problems  of  gas  dynamics,  plasma  physics  and  astrophysics  [1.9- 
10].  In  particular  it  gave  reliable  results  at  solution  of  a  problem  of  nitrogen  plasma  jet  interaction 
with  oncoming  flow  [1.19]. 
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Chapter  2.  Impact  of  plasma  from  cylindrical  plasma  generator  on  counter  and  cross  flows 

2.1.  Approach  to  solution  of  a  problem  of  impact  of  plasma  from  cylindrical  plasma 
generator  on  flow. 

Before  clarifying  which  parameters  of  plasma  generator  with  the  divergent  nozzle  are 
preferable  we  made  modeling  of  impact  of  plasma  from  cylindrical  plasma  generator  on  counter  and 
cross  flows  for  known  experimental  parameters.  This  allowed  to  test  the  opportunity  of  application 
of  developed  numerical  method  to  problems  of  plasma  interaction  with  counter  and  cross  flows. 

We  also  considered  the  interaction  of  flat  jet  with  the  supersonic  cross  flow  as  in  [2.1].  In 
Fig.2.1  one  can  see  the  principle  scheme  of  such  interaction  of  the  jet  with  the  crossflow.  The 
crossflow  consisting  of  a  flammable  mixture  interacts  with  plasma  jet  at  some  angle  to  a  flow  (0  =  0 
-180°). 

As  in  our  simulation  in  Ref  1,  the  gasdynamic  equations  are  coupled  with  the  2-D  Maxwell 
equations  and  the  plasma  dynamics  equations,  accounting  for  ohmic  heating  and  radiative  plasma 
cooling.  All  conditions  that  are  modeled  correspond  to  the  experiments  in  Refs.  2,3.  Specific 
parameters  are  discharge  current  1=  45  A,  Nitrogen  mass  flow  mf  =  0.7  g/s,  power  input  P=7.5  kW 
and  a  freestream  Mach  number,  M=2.  Appropriate  thermodynamic,  optical  and  transport  properties 
of  the  nitrogen  plasma  are  used  in  the  Computations.  A  Lagrange  method  based  on  a  completely 
conservative  implicit  difference  scheme  with  an  adaptive  triangular 
unstructured  grid  described  in  Chapter  1  is  applied.  This  was  effective  in  modeling  the  coupling  of 
the  internal  and  external  flows  inside  and  outside  the  counterflow 
jet  and  its  nozzle  [2.1].  It  effectively  simulated  special  flow 
features  that  complicated  the  modeling.  Some  of  these  were  1) 
greatly  disparate  length  scales  that  were  of  the  order  of  1  mm 
inside  the  plasma  generator  and  100  mm  in  the  external  flow 
stream-tube  capture  cross  section,  2)  mixed  subsonic  and 
supersonic  flow  regions,  3)  large  gradients  of  plasma 
conductance  (over  3  decimal  orders  of  magnitude  for  1  mm),  and 
4)  other  large  gradients  requiring  fine  grids  to  resolve. 

In  what  follows,  this  approach  is  applied  in  this  investigation  of 
interaction  of  plasma  jets  with  crossflows. 


i-t 

a> 

a> 
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Fig.2. 1  Plasma  jet  from  a  slot  created  by 
a  row  of  plasma  generators. 
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2.  2.  Plasma  generators  injecting  plasma  jet  to  a  crossflow 

Approach 

The  model  is  based  on  the  following  system  of  radiating  plasma  electrodynamics  and  gasdynamic 

equations  in  more  details  it  is  represented  in  Chapter  1 . 

dpidt  +  div(pv)  =  0,  (2.1) 

p  dv/dt  =  -  grad  p, 

p  deldt  =  -  p  div  v  +  j^/a  -  q, 

rot  E  =  - 

rot  B  =  poj, 

div  B  =  0, 

j  =  aE, 

where  p,  p,  v,  s  are  respectively  the  plasma  (gas)  density,  pressure,  velocity,  and  internal  energy;  a 
is  the  electric  conductance  of  the  gas  and  plasma,  q  is  the  specific  power  of  plasma  thermal 
radiation,  E  is  the  electric  field  intensity,  B  is  the  magnetic  field  induction,  j  is  the  electric  current 
density,  and  po=47r  lO“^  H/m.  An  unsteady  formulation  was  applied  in  cylindrical  coordinates. 
Accordingly,  v  =  (vr,  Vjp,  Vz),  B  =  (0,  B(p,  0),  j  =  (jr,  0,  jz),  E  =  (Er,  0,  Ez),  d/d(p  =  0;  the  subscripts  r, 
(p,  z  correspond  to  the  radial  r,  azimuth  (p,  and  axial  coordinate  z.  The  plasma  equation  of  state  can 
be  expressed  as 

p  =  a phsT/nimoi,  s=p/[p{r-'^)],  (2.2) 

where  yis  the  adiabatic  exponent,  ks  is  the  Boltzmann  constant,  T  is  the  plasma  temperature,  nimoi  is 
the  initial  gas  molecular  mass,  and  the  factor  a  accounts  for  the  free  species  number  changes 
resulting  from  chemical  reactions  (primarily  dissociation  and  ionization). 

The  functions  y  =  y  (p,T),  a  =  a(p,T),  a  =  a  (p,T)  and  q  =  q(p,T)  for  nitrogen  plasmas  were 
computed  estimated  initially  from  data  in  Refs.  4-6.  Calculations  first  assumed  a  =  1.1  and  y  =1.4, 
but  changed  with  subsequent  parametric  variations  of  the  temperature,  T. 

The  computational  region  (see  Fig.  T.2)  includes  both  the  regions  inside  and  outside  the 
plasma  generator.  No-slip  boundary  conditions  on  the  walls  were  assumed  for  both  plasma  and  gas 
and  d/dr  =  0  was  applied  on  the  axis  of  symmetry. 

The  gas  input  azimuthal  and  radial  velocities  v^,  Vro  were  determined  from  the  gas  density 
po  of  the  given  gas  flow  rate  m'  and  areas  of  the  inflow  tangent  orifices  Fj,  (see  Fig.  2),  and  of  the 
circular  input  slot  v^,  =  m'/{pH)  F^,  Vro=  m'/{po  Fg). 
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3  p  /d  r  =  0  , 


Fig.  2.2  Flow  legions;  1  -  gas  injection  slot,  2  -  arc,  3  -  plasma  generator's  chamber,  4  -  outlet  nozzle, 
5  -  external  region. 

The  input  gas  temperatme  To  =  lb?"  K.  Conditions  p  =  pf,  Vz^-Vrf  ,  T  =  Tf  on  the  right-hand 
boundaiy  defined  the  supersonic  ambient  flow.  Conditions  on  the  left-hand  boimdaiy  and  on  the 
peripheral  boundary  provided  approach  to  the  ambient  fteestream  values. 

The  entire  left-hand  wall  inside  the  plasma  generator  had  an  electric  potential  in  the  slot 
where  the  gas  is  injected,  the  condition  on  the  magnetic  field  was  defined  fiom  the  total  discharge 
cuiient  1.  Accordingly,  Bo  =  poI/(2Trro),  where  tq  is  the  distance  between  the  axis  and  the  wall  of  the 
jet  channel.  All  the  remaining  walls  were  at  the  electiic  potential  qyi.To  expedite  the  computations, 
a  small  background  conductance  a  <  10”^  (Tmax  was  assumed  to  exist  in  the  cold  gas;  here,  <7„.av  is 
the  highest  value  of  plasma  conductance  in  the  electiic  cmrent  chaimel.  The  computed  distiibution 
of  electric  ciurent  in  the  rest  of  the  channel  proved  to  be  practically  independent  of  the  fluctuations 
of  the  plasma  conditions.  This  is  caused  by  the  weak  dependence  of  the  nitrogen  plasma 
conductance  on  plasma  pressure  and  temperatme  in  the  parametric  region  of  interest.  Accordingly, 
it  was  possible  to  “freeze”  this  distribution  and  to  avoid  its  expensive  recalculation  during  the 
computations.  All  the  data  correspond  to  the  conditions  in  Refs  4,6.  There,  the  number  of  gas 
inflow  orifices  was  fom  and  their  diameters  were  1  mm. 

As  shown  in  Ref.  2,  the  distribution  of  temperature,  pressure  and  density  varies  over  the 
nozzle’s  exit.  This  showed  the  main  difference  between  our  computations  and  those  of  Refs.  4  and 
6,  where  these  disftibutions  were  assumed  constant  fiuictions.  The  distinction  between  these 
calculations  and  Refs.  2.3  and  2.5  is  connected  with  internal  heat  transfer  processes  from  the  source 
of  energy  release  inside  the  plasma  generator. 
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Fig.  2.3.  Computational  grid  inside  the  plasma  generator. 

In  Fig.2.3-7  are  represented  distributions  of  parameters  inside  and  outside  of  the  plasma  generator, 
when  injection  of  plasma  takes  place  at  the  angle  0=180°.  In  Fig.2.3-4  one  can  see  the 
computational  grid  inside  and  outside  of  the  plasma  generator. 


z 


Fig.  2.4.  Computational  grid  outside  the  plasma  generator. 

hi  Fig.2.5  one  can  see  computation  results  results  of  pressure  distribution  inside  and  outside  the 
plasma  generator  (in  imits  10^  Pa)  for  time  moments  t=2.93  ms  (a),  4.87  ms  (b)  h  1.57  ms  (c). 


Fig.  2.5a 
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Fig.  2.5b 


Fig.  2.5c 


Fig.2.5.  Pressure  distributions  (in  10^  Pa)  outside  and  inside  of  the  plasma 
generator  for  time  moments  t=2.93  ms  (a),  4.87  ms  (b)  h  1.57  ms  (e) 


In  the  Fig.2.6  one  ean  see  distributions  of  temperature  (in  units  T/Tq)  outside  and  inside  of  the 
plasma  generator  for  the  same  time  moments. 
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Fig.2.6  a 
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Fig.2.6  b 


Fig.2.6  c 


Fig.2.6.  Temperatuie  distributions  (in  T/To)  outside  and  inside  of  the  plasma 
generator  for  time  moments  ?=2.93  ms  (a),  4.87  ms  (b)  h  1.57  ms  (c) 


In  Fig.2.7.  one  can  see  distributions  of  the  absolute  value  of  velocity  outside  and  inside  of  the 
plasma  generator  for  the  same  time  moments  (in  units  10^  m/s). 


Fig.  2.7  a 
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Fig.  2.7  b 
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Fig.  2.7c 
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Fig.2.6.  Absolute  value  of  veloeity  distributions  (in  10  m/s)  outside  and 
inside  of  the  plasma  generator  for  time  moments  t=2.93  ms  (a),  4.87  ms  (b)  h 
1.57  ms  (e) 

Condueted  eomputations  allow  to  determine  distributions  of  main  parameters  at  the  plasma 
generator  exit  (cut)  necessary  for  modeling  of  interaction  of  plasma  jet  with  a  flow  in  case  of  the 
plasma  jet  incidence  to  the  flow  0  =  0  -180°. 

Fig.2.5-7  show  that  the  plasma  jet  significantly  effects  the  flow.  Formation  of  local  circular 
regions  with  higher  temperature  and  lower  pressure  take  place  in  the  flow.  Favorable  conditions  for 
ignition  of  a  fuel  mixture  are  created  here.  However  this  geometry  of  interaction  is  the  model  one  ( 
the  plasma  is  directed  along  the  wall  of  a  combustor)  because  the  interaction  of  the  plasma  with  the 
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wall  of  the  combustor  would  sooner  destroy  the  latter  due  to  high  temperature,  so  it  can  be 
considered  as  the  necessary  step  for  calculation  of  plasma  flow  interaction  at  other  angles. 

One  can  see  from  the  figures  that  creation  of  vortices  and  stagnation  zones  takes  place  in  the  cavity 
of  the  «barrel»  type  inside  the  plasma  generator,  it  according  to  [2.6]  can  lower  the  temperature  at 
the  exit  from  the  plasma  generator,  and  respectively  to  decrease  efficiency  of  the  plasma  jet 
interaction  with  the  flow.  So  in  computations  for  the  aims  of  the  Project  we  have  conducted 
analyses  of  the  conical  channel  of  plasma  generator  design  (see  below). 

To  investigate  the  influence  of  the  radial  distribution  of  the  flow  variables  over  the  nozzle  exit,  we 
made  two  series  of  computations.  The  first  used  constant  flow  conditions  at  the  nozzle’s  exit  and 
the  second  allowed  for  variations.  For  this  study,  the  plasma  generator’s  radius  r=W  m,  pressure 
/)=0.235  atm,  temperature  r=167°  K,  density  p=0.5  kg/m  ,  Mach  number,  M=1  was  assumed.  In 
addition  the  jet  is  directed  against  the  flow  (0=180°)  (counterflow  case).  The  characteristic  plasma 
flow  time  in  the  calculations  is  about  10'"^  s.  To  resolve  the  large  gradients  within  the  nozzle  and 
near  the  plasma  generator,  the  total  number  of  calculated  points  was  about  20,000^. 

Calculations  for  plasma  jets  injected  at  angles  0  =  135°  and  90°  with  respect  to  the  free  stream 
crossflow  have  been  performed.  As  contrasted  to  the  internal  nozzle  flow,  the  total  number  of  grid 
points  for  the  external  flow  field  was  only  about  8,000.  However,  to  resolve  the  strong  gradients, 
the  grid  variability  parameter  was  100.  The  jet  initial  profile  was  considered  to  be  “top  haf’ 
(constant  over  the  jet  exit  region  and  zero  off  it),  modeling  a  row  of  equally  spaced  plasma 
generators  shown  in  Fig.  2.1.  This  formulation  is  intended  to  clarify  characteristics  of  the  jet 
interaction  with  the  crossflow,  spatial  distribution  of  the  flow  variables,  distribution  and  special 
features  of  the  jet  influence  on  the  oncoming  flow  including  formation  of  shocks,  dependence  on 
outflow  angles,  etc.  To  expedite  our  initial  modeling  approach,  we  ignore  viscosity  and  heat 
conduction  effects  that  are  generated  on  length  scales  of  the  order  of  0.01  mm.  These  are 
infinitesimal  on  our  computational  domain.  We  neglect  the  global  effect  of  these  processes  in  this 
analysis  as  well  as  combustion  and  turbulence  phenomena. 

The  free  Lagrange  method  based  on  a  completely  conservative  implicit  difference  scheme  as 
well  as  adaptive,  automatically  unstructured  triangular  girding  was  applied  in  the  numerics.  The 
characteristic  size  of  the  cells  near  the  exit  nozzle  (where  the  plasma  jet  is  situated)  was  about  0.1 
mm.  This  was  found  to  be  appropriate  to  simulate  large  and  medium-size  vortices  and  corresponding 
heat,  momentum,  and  mass  transfer  processes 
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2.3.  Plasma  properties 

In  order  to  obtain  initial  results  on  plasma  eomposition  we  started  from  Refs.  7  and  8  for  nitrogen 
equilibrium  thermal  plasmas.  These  showed  that  N2  molecules,  N  atoms  and  atomic  ions  play 
key  roles  in  the  temperature  range  of  interest,  which  is  5000-10000°K.  These  species  participated  in 
the  following  dissociation  and  ionization  reactions: 

N+N^N^, 

Concentrations  of  atomic  ions  at  these  temperatures  are  much  greater  than  those  of  molecular 
N2^  ions  that  appear  in  the  reaction 
e  +  Nt  -^N^. 

The  corresponding  equations  for  these  reactions  are  the  standard  ones  Ref.  9: 

N^/N,=Pl(T),  (2.3) 

(N*)- /N  =  F2(T),  (2.4) 

and  also  the  equation  of  mass  balance 

p/mj^=N  +  2N2 , 

where  p  is  nitrogen  density,  mN-weight  of  the  nitrogen  atoms.  Thus, 

FI  =  3.34-10^^  •(!- 

exp(-3400  /  T))  ■  exp(-l  1 3265  /  T),  cm  ' 

F2  =  1 .22  •  1  O'  •  r* '  •  exp(- 1 68677  /  T),  cm'' . 

These  equations  were  used  for  our  plasma  modeling. 

2.4.  Discussion  of  plasma  jet  results 

Figure  2.8  shows  our  experimentally  obtained  flow  patterns  for  a  hydrocarbon  plasma  jet  interacting 
with  an  air  crossflow  for  ^=90°.  The  velocity  and  temperature  of  the  plasma  jet  were  Vj«850  m/s, 
T«6000°K  respectively.  The  Mach  number  and  temperature  of  the  crossflow  were  M=2,  T=167°K. 
These  characteristics  are  close  to  those  of  experiments  [2.3,  2.5]  whose  conditions  we  used  in  our 
computations.  Figure  2.8  shows  that  the  plasma  jet  is  strongly  deflected  by  the  crossflow.  Moreover, 
two  regions  are  formed  that  have  distinctly  different  temperatures.  One  (red-moderately  hot)  part  is 
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long  and  the  other  (blue  and  white-very  hot)  part  is  thin.  This  structure  is  consistent  with  our 
computations  that  are  shown  in  Figs.  T9  to  be  discussed  in  what  follows. 


Fig.  2.  8.  Time-averaged  photo  exposure  from  our  experiment  of 
plasma  jet  interaction  with  the  crossflow  at  0=  90°.  Velocity  of  the 
plasma  jet  Vj »  850  m/s,  T»6000  K.  Mach  number  of  the  crossflow 
M=2,T=167°  K. 


Figure  2.9  corresponds  to  a  cross  section  of  the  jet  flow 
injected  at  ^  =  135°  into  nitrogen.  In  Fig.  2.10,  isotachs  are 
shown  for  6  =  90°  injection  into  nitrogen.  Characteristics  of 
the  crossflow  are  the  same  as  previously  considered. 


It  is  apparent  that  for  the  conditions  investigated  in  this  paper,  the  computations  for  injection 
angles  of  135°  and  90°  give  similar  results.  This  suggests  that  at  least  for  our  parameter  space  the 
effect  of  injection  angle  may  be  weak.  This  observation  may  apply  to  combustion  processes  in  more 
general  parametric  ranges.  If  this  conjecture  is  correct,  it  can  be  exploited  to  reduce  the  number  of 
computations  in  which  reactive  species  are  assessed.  The  results  obtained  for  these  cases  are  in 
qualitative  agreement  with  the  experimental  picture  of  the  flow  in  Fig.  2.8.  In  contrast,  results  for 
uniform  conditions  at  the  nozzle  exit  differ  considerably  from  experimental  and  non-uniform  ones. 

The  computations  show  that  during  time  intervals  of  the  order  of  (5-10)xl0''^  s,  the  jet 
expands  and  accelerates  from  the  nozzle  face.  During  this  period,  the  Mach  number  reaches  M~3 
even  in  the  hot  part  of  jet,  where  T~4000°  K.  In  addition,  the  jet  screens  the  flow  and  a  bow  shock  is 
formed  ahead  of  it.  The  oncoming  flow  is  deflected  by  the  shock  to  flow  along  the  jet.  These 
processes  result  in  velocity  differences  between  the  jet  and  surrounding  flow  as  much  as  3000  m/s. 

Mixing  of  the  fluid  from  the  jet  occurs  with  the  region  it  shadows  behind  itself  as  well  as 
that  ahead  of  itself  Low  values  of  pressure,  temperature  and  density  in  the  cavity  behind  the  jet. 
Some  vortical  structures  are  observed  on  the  mixing  layers,  possibly  due  to  computational  viscosity 
in  this  nominally  inviscid  formulation.  These  are  consistent  with  those  observed  in  the  experiments 
discussed  in  Refs.  2.10,  2.12  and  2.13.  A  cavity  is  formed  behind  the  jet  and  its  size  increases  with 
time.  In  a  real  fluid,  this  cavity  would  represent  a  separated  flow  region  in  which  the  vorticity  is 
driven  by  the  jet  shear  layer  entrainment.  References  2.10,  2.12  and  2.13  show  coherent  large  scale 
vortical  structures  from  such  shear  layers.  Expansion  of  the  cavity  is  accompanied  by  creation  of  a 
bow  shock. 


The  temperature  plots  in  Fig.2.9d  show  patches  of  heated  gas.  They  resemble  experimentally 
observed  elongated  structures  on  the  surface  of  a  hydrogen  jet  injected  into  a  supersonic  crossflow 
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discussed  in  Ref.  2.12.  More  detailed  calculation  of  the  shear  layer  between  jet  and  oncoming  flow 
shows  Rayleigh-Taylor  instabilities  Ref  2.1 1.  The  size  of  area  where  T>900°  K  is  much  larger  than 
that  of  the  dissociated  gas  as  shown  in  Figs.  2.9-11.  Mixing  associated  with  this  instability  could 
facilitate  combustion  processes  inside  a  combustor. 

The  density  variations  in  Fig.  2.9c  show  that  the  jet  bow  shock  travels  with  a  velocity 
V~  1300  m/s  at  the  time  t=0.5  -10'"^  s.  This  is  associated  with  a  shock  Mach  number  M~3. 

Figure  2.11  shows  equilibrium  concentration  contours  of  atoms  N  and  ions  obtained  by 
formulas  (2.3)-(2.4).  It  is  evident  that  the  concentrations  of  atoms  are  high.  They  are  the  same  order 
as  those  of  the  molecules  for  which  N~10**  cm'^  in  the  jet.  Concentrations  of  ions  are  high  -lO'^ 
cm'  .  The  region  that  they  occupy  is  much  smaller  than  that  of  the  atoms.  It  also  appears  that  the 
zone  of  atom  production  is  much  smaller  than  that  of  the  gas  heated  region.  For  air  however,  it  will 
be  much  larger,  since  the  dissociation  energy  of  oxygen  molecules  is  about  one  half  of  that  of  the 
nitrogen  molecules.  This  can  enhance  combustion  and  propulsion  efficiency. 


Pressure  (1  atm),  time=0.57*10'^  s. 
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Fig.  2.9a  Isomachs  (0=135°) 


Fig.  2.9b  Isobars  (0=45°) 
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Density  (10  kg/m^),  time=0.57*10'*  s. 
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Fig.  2.9c  Density  contours  (0=135°) 


Temperature  (100  K),  time=0. 57*10"^  s. 
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Fig.  2.9d  Isotherms  (0=45°) 
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Fig.  2.  9e  Isomachs  (0=135°) 
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Speed  (100  m/s),  time=0.57*10'*  s. 
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Fig.2.10  Plasma  jet  interaction  with  the  crossflow  at  0=90°:  T=3000  K,  M=l,  P=3  atm.  The 
crossflow  parameters:  M=2,  P=0.25  atm,  temperature  T  =167  K.  Flow  speed  contours. 
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Fig.  2.11a  Spatial  distribution  of  atoms  in  nitrogen  (N,  10^^  cm’^)  (0=135°) 
excited  by  nitrogen  plasma  jet. 
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Fig.  2.12.1)  Spatial  distribution  of  atomic  ions  in  nitrogen  (N^,  10^^  cm'^)  (6=135°) 
excited  by  nitrogen  plasma  jet. 

2.7.Coiiclusions  to  Chapter  2. 

This  chapter  presents  a  theoretical-computational  analysis  of  plasma  jet  interaction  with  a 
crossflow.  It  emphasizes  special  features  of  the  coupling  of  the  plasma  chemistry  and  gasdynamics 
of  importance  to  ignition  and  combustion  in  hypersonic  propulsion  systems.  This  analysis  is  a  first 
step  to  assess  the  effectiveness  of  plasma-jet  injection  for  augmentation  of  ignition-combustion 
processes. 

Equations  for  plasma  jets  in  a  crossflow  including  transport,  chemical  and  ionization  processes 
were  applied.  The  equations  derived  were  numerically  treated  by  an  implicit  fiee-Lagrange  method 
using  and  adaptive  unstructured  grid. 

Om'  computations  show  a  complex  and  non-uniform  structme  of  the  flow  outside  the  plasma 
generator.  High  atom  concentrations  inside  the  plasma  jet  are  apparent.  This  fact  differentiates 
plasma  jets  from  hot  gas  jets  and  can  lead  to  combustion  and  propulsion  enhancement  over  that 
possible  with  the  latter.  Oiu'  results  are  in  qualitative  agreement  with  our  own  experiments  and 
those  of  others.  It  is  hoped  that  these  results  will  stimulate  new  studies  that  clar  ify  the  relative  roles 
of  plasma  and  gasdynamic  processes  in  these  flows. 
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Chapter  3.  Computations  of  gasdynamic  parameters  of  plasma  generator  with  divergent 

channel 


3.1.  Formulation  of  a  problem. 

Starting  from  the  interest  of  analysis  of  application  possibility  of  plasma  generator  with 
divergent  channel  for  influencing  on  supersonic  gas  gaseous  mixure  flows  we  carried  out 
computations  for  the  plasma  generator  represented  in  experimental  section  of  the  Report  (see 
below). 

From  a  gasdynamic  point  of  view  such  a  plasma  generator  is  a  generator  with  a  conical 
nozzle  [3.1-2],  The  detailed  gasdynamic  analysis  of  conical  nozzle  is  presented  in  [3.3],  according 
to  it  a  flow  in  such  a  nozzle  represents  a  complicated  structure.  In  vicinity  of  a  conjunction  point 
between  the  radius  and  conical  parts  of  a  nozzle  can  take  place  a  flow  stagnation  at  some 
conditions.  A  compression  wave  origins  a  shock  wave.  It  reflects  many  times  from  a  symmetry  axis 
and  of  the  nozzle  wall.  Its  intensity  increases  at  approaching  to  the  axis  of  symmetry,  and  farther 
down  the  flow  it  decreases  because  of  an  influence  of  rarefaction.  In  majority  of  cases  an  intensity 
of  a  wave  is  not  great.  Existence  of  a  shock  leads  to  non-homogeneity  in  distribution  of  a  flow 
parameters  over  cross  sections  of  a  flow. 

Parameters  of  the  plasma  generator  under  the  designing  are  the  follows:  Minimal  nozzle 
diameter  3  mm.  Full  angle  12°,  half  angle  6°.  Nozzle  length  along  the  axis  30  mm.  Cathode 
diameter  4.5  mm.  Half  angle  of  the  conical  cathode  75°.  The  distance  from  the  cathode  to  the 
nozzle  1.7  mm.  An  area  of  the  pass  in  channels  equals  to  an  area  in  the  critics  of  the  nozzle. 
Nominal  power  10  kW,  gases  nitrogen  or  argon,  energy  input  is  distributed  at  a  distance  of  5  mm  in 
front  of  the  cathode. 

In  Fig  3.1  one  can  see  the  computation  grid  for  such  a  plasma  generator  used  in  our 
calculations. 

The  calculation  data  is  the  follows.  Minimal  nozzle  diameter  3  mm.  Full  angle  12°,  half 
angle  6°.  Nozzle  length  along  the  axis  30  mm.  Cathode  diameter  4.5  mm.  Half  angle  of  the  conical 
cathode  75°.  The  distance  from  the  cathode  to  the  nozzle  1.7  mm.  An  area  of  the  pass  in  channels 
equals  to  an  area  in  the  critics  of  the  nozzle.  Nominal  power  10  kW,  gases  nitrogen  or  argon,  energy 
input  is  distributed  at  a  distance  of  5  mm  in  front  of  the  cathode. 

The  calculated  plasma  generator  with  nitrogen  as  plasma  forming  gas  variant  has  the 
following  parameters:  gas  consumption  -  0.7  g/s.  Electric  power  of  the  source  9  kW.  Ratio  of 
translational  and  rotational  (at  inlet)  velocities  equals  to  tg(6°).  Plasma  generator  injects  into  a 
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normal  medium  at  pressure  1  atm  and  temperature  300  K.  At  the  starting  moment  in  the  cross 
section  of  the  inlet  conditions  are  normal:  P=1  atm  and  T=300  K. 


Fig.3.1.  Computational  grid  of  internal  and  external  regions  of  the  plasma  generator  being  under  the 
design. 


3.2  Initial  calculations.  Initially  calculations  have  been  conducted  only  for  inflowing  of  a 
gas  in  the  plasma  generator  in  an  absence  of  energy  input  to  a  discharge,  conditions  of  so  called 
“cold”  plasma  generator.  In  Fig.  1-4  are  presented  the  data  on  pressure  (p  is  in  atm),  temperature  (T 
in  100  K),  Mach  number  (vmah=  I  v  I  / 1  a  (local)  I )  and  z  component  of  a  velocity  (vz  in  100  m/s) 
distributions  in  a  «cold»  plasma  generator.  Indicated  time  corresponds  to  stationary  parameters  to 
which  the  flow  has  approached. 

Calculations  show  that  some  increase  of  pressure  in  the  plasma  generator  conical  part  takes 
place  near  walls,  where  the  gas  heating  takes  place.  Formation  of  the  vortex  area  takes  place  near 
the  channel  center  in  the  vicinity  of  plasma  generator  outlet.  Distributions  of  vz  velocity 
component  and  of  Mach  number  outside  the  plasma  generator  show  appearance  of  complex  vortex 
structures  at  the  exit  from  the  plasma  generator. 
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Fig.  3.2,  3.  Distributions  of  pressure  and  temperature  in  a  «cold»  plasma  generator. 

Pressure  (p  in  atm),  temperature  (T  in  100  K),  z,  r  -  in  cm. 
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Fig.  3.  4,5.  Distributions  of  Mach  number  module  and  of  z  component  of  a  velocity  in  a  “cold”  plasma  generator. 

vz  in  100  m/s  ,  M=vmach=  |  v  I  /  |  a  (local)  I ,  z,  r  -  cm. 
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3.3  Comparison  with  available  experiments. 

In  calculations  we  located  the  zone  of  energy  input  in  the  area  near  the  eleetrode  (at  a  distanee 
~L/3,  where  L  is  a  length  of  the  plasma  generator)  it  produeed  the  energy  input  with  the  power 
equal  to  those  of  the  plasma  generator  as  it  was  done  for  acceleration  of  caleulations  in. 

In  Fig.  3.6.a-d  results  of  plasma  generator  work  dynamics  in  time  are  presented,  so  ealled 
“hot”  plasma  generator.  They  show  that  the  stationary  state  of  the  plasma  generator  is  realized  not 
immediately  but  in  the  result  of  several  intermediate  states  leading  to  reeonstruction  of  a  flow.  After 
that  the  state  with  high  temperature  in  the  eenter  of  a  ehannel  is  reached. 

Obtained  results  are  in  qualitative  agreement  with  experimental  ones  obtained  for  analogous 
plasma  generator  [3.4-5]  with  the  eonieal  diverging  plasma  nozzle  and  parameters  (diameter  5-6 
mm  of  the  anode,  diverging  angle  6°,  eonsumption  1.5  g/s  and  power  50-100  kW),  which  by 
minimum  power  is  by  ~  5  times  greater  than  our  plasma  generator  and  by  nitrogen  eonsumption  is 
greater  by  ~  2  times. 

The  temperature  in  T  [3.5]  in  initial  seetion  ~  L  /3  changed  over  radius  from  0.3  to31  kK, 
and  in  our  eomputations  from  0.3  to  ~45  kK  and  had  elose  distribution  over  radius  to  those  of  [3.5] 
(see  Fig.  3.6c.  Higher  values  in  our  computations  in  comparison  with  [3.5]  could  be  eonditions  by 
negleet  of  thermal  eonductivity  processes  in  computations,  their  aeeount  should  lead  to 
redistribution  of  temperature  with  deereasing  of  its  value  in  a  center  of  a  ehannel.  Values  of  mean 
mass  density  of  the  same  order  of  magnitude  (5-6  g/cm  ),  we  also  observed  analogous  radial 
distribution  of  gas  density. 

Electron  concentration  was  about  Ne»8-10  cm'  .  Electron  concentration  in  [3.5]  was 
proposed  a  meehanism  of  additional  photoionization  due  to  VUV  radiation  of  the  plasma  for  a 
description  of  such  a  high  values  of  eleetron  eoncentrations.  High  values  of  eleetron  eoneentrations 
in  our  computations  are  conditioned  by  the  high  gas  temperature. 

Clarification  of  real  reasons  of  this  difference  requires  additional  experimental  and 
theoretical  investigations  with  the  definite  eonditions  in  the  arc  produced  by  our  plasma  generator. 

Computation  results  on  the  atom  eoncentrations  give  large  values  of  the  order  N«5-7  -lO'^ 
em'^,  this  shows  the  prospeets  of  this  plasma  generator  application  for  the  problems  of  a  fuel 
ignition. 

Condueted  analysis  and  eomparison  with  experiment  show  principally  correct  modeling  of 
processes  in  the  plasma  generator  of  the  given  type. 
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Fig.  3.6  a  Pressure,  temperature,  vz  and  vmach  distributions  in  a  “hot”  plasma 
generator  at  IS  IO"^  s. 

Pressure  (p  ,  atm);  temperature(T  ,100  K);  vz  ,  100  m/s  ,  vmach=  |  v  I  / 1  a  (local)  I  z,  r  -cm. 
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Fig.  3.6b  Pressure,  temperature,  vz  and  vmach  distributions  in  a  “hot”  plasma  generator  at 
19.5-10''‘  s.  Pressure  (p  ,  atm);  temperature(T  ,100  K);  vz  ,  100  m/s  ,  vmach=  I  v  I  / 1  a  (local)  I  z,  r  - 
cm. 
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Fig.  3.6c.  Pressure,  temperature,  vz  and  vmach  distributions  in  a  “hot”  plasma  generator  at 
28.3-10''*  s.  Pressure  (p  ,  atm);  temperature(T  ,100  K);  vz  ,  100  m/s  ,  vmach=  |  v  I  / 1  a  (loeal)  I  z,  r  - 
cm. 
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Eig.  3.6d  distributions  of  electron  Ne  and  atom  N  concentrations  in  a  “hot”  plasma  generator 
at  28.3*10^  s.  Ne  is  in  10^^  cin'^ ,  N  is  in  lO’®  cm'^,  z,  r  -cm. 


3.4.  Investigations  of  interaction  of  a  jet  with  the  oncoming  flow 

Dming  the  next  stage  of  investigations  we  plan  to  study  the  action  of  the  nitrogen  jet  from 
the  conical  nozzle  plasma  generator  on  the  oncoming  airflow.  As  showed  om  previous  calculations 
it  is  necessary  to  know  the  distribution  of  its  main  parameters  at  its  outflowing  at  6  =180°  into  a 
gas.  So  we  have  imdeitaken  such  computations  of  interaction  of  the  plasma  generator  jet  with 
nifrogen  flow  at  6=180°. 
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We  have  undertaken  our  eomputations  for  eounter  flow  with  a  Mach  number  M=2.  In  Fig. 
3.7  a-f  one  can  see  computation  data  of  interaction  of  a  such  plasma  jet  with  nitrogen  flow  not  long 
before  its  coming  to  a  quasi-stationary  state. 


Fig.  3.7  a,b  Pressure  and  temperature  distributions  in  a  plasma  generator  jet  interacting  with 
the  oncoming  flow  at  the  incidence  angle  0=180°.  Nitrogen,  M=2,  Poo=  1  atm,  Too=  300  K,  z,  r  -cm. 
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Vz,  100  m/s,  time=  1.2*1  O' s 
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Fig.  3.7  c,d  Distributions  of  a  Mach  number  and  Vz  velocity  in  a  plasma  generator  jet 
interacting  with  the  oncomhig  flow  at  the  incidence  angle  0=180°.  Nitrogen,  M=2,  Pqo=  1  atm, 


Too=  300  K,  z,  r  —cm. 


N„  ■  cm  ,  time=1 .2*1 0  s 


Fig.  3  .7  e,f  Dishibutions  of  election  Ne  and  atomN  concentrations  in  a  plasma  generator 
jet  mteracting  with  the  oncoming  flow  at  the  incidence  angle  6=180°.  Nitiogen,  M=2,  ?<»=  1  atm, 
Tod=  300  K,  z,  r  —cm. 
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Computations  have  been  carried  out  in  a  following  way.  Start  of  computations  with  a 
source  of  a  disturbance  (counterflow)  took  place  at  time  moment  2.08  s  where  gasdynamic  data 
obtained  on  plasma  jet  injection  to  a  motionless  gas  (see  above)  have  been  used.  The 
counterflow  to  the  stream  was  given  at  z>10  cm.  This  lead  to  a  decay  of  the  discontinuity.  A 
shock  goes  forward  with  increased  pressure  (more  than  3)  and  Mach  number  (less  than  0.5).  A 
back  wave  in  a  decay  of  the  discontinuity  is  also  a  shock  one.  It  moves  slowly  along  a  flow.  On  a 
left  outletting  boundary  we  have  zero  boundary  condition  of  the  outletting  velocity,  so  on  this 
boundary  a  shock  also  has  appeared.  In  the  next  moment  we  changed  the  boundary  condition  on 
the  outletting  boundary.  We  put  there  the  outletting  velocity  generating  the  rarefaction  wave  up 
the  flow,  it  was  accepted  to  be  equal  to  a  velocity  of  the  oncoming  flow  in  order  to  eliminate  the 
shocks  of  density  increase  in  the  oncoming  flow. 

In  connection  with  these  results  we  now  have  a  computational  information  about  the  plasma 
generators  operation  in  the  subsonic  flow  of  the  increased  pressure.  Its  peculiarities  are:  the 
oncoming  flow  locks  the  plasma  generator,  a  gas  goes  out  only  in  a  thin  layer  along  a  surface  of 
a  nozzle.  The  pressure  in  the  plasma  generator  is  increased,  the  radiation  increases  also 
(proportionally  to  square  of  the  pressure)  and  the  summed  effect  of  the  gas  heating-expansion 
becomes  smaller. 

Now  we  conduct  continuous  computations  of  penetration  of  a  jet  through  the  oncoming  flow 
until  it  comes  to  a  stationary  state.  Obtained  corresponding  values  will  be  used  in  solution  of  a 
problem  of  air  stream  activation. 

Our  computations  have  shown  that  the  distribution  of  flow  parameters  at  the  outlet  of  the  plasma 
generator  has  highly  non  homogeneous  character.  For  smoothing  of  these  distributions  we  have 
to  account  effects  of  thermal  conductivity  and  viscosity  in  the  plasma  generator  channe  data  on 
which  is  represented  in  [3.6-7]. 

3.5.  Viscosity,  thermal  conductivity  and  different  angles  of  swirling  inside  the  plasma 

generator  influence  on  plasma  parameter  distributions 

In  this  part  of  the  report  we  present  computation  results  of  the  plasma  jet  injection 
towards  the  flow.  The  purpose  of  computations  was  to  clarify  the  viscosity  and  thermal 
conductivity  influence  and  the  different  angles  of  swirling  inside  the  plasma  generator  on  the 
plasma  parameter  distributions  over  the  out  letting  nozzle  cross  section.  These  parameter 
distributions  [3.2]  later  are  applied  for  modeling  of  plasma  jet  interaction  with  a  cross  flow  at  the 
angle  (18O°>0>O°). 
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The  computation  jet  variants  corresponded  to  the  jet  injection  to  the  oncoming  supersonic 
flow  with  the  parameters:  the  incident  flow  velocity  is  M=2,  the  pressure  0.2  atm.  and  the 
temperature  is  250°  K. 

As  earlier  the  mass  flow  rate  was  given  at  the  inflow  to  the  plasma  generator,  the  pressure 
is  extrapolated  inside  the  plasma  generator  to  the  inflow  boundary,  density  and  the  temperature 
are  determined  from  the  condition  of  the  isentropic  gas  transition  from  the  normal  state  to  the 
state  with  the  determined  pressure  ,  the  velocity  of  the  inflow  is  determined  by  the  given  mass 
flow  rate. 

The  figures  in  all  the  computation  variants  are  represented  for  the  time  moment  when  the 
plasma  jet  head  reached  approximately  15  cm. 

Computation  variants  are  the  following: 

B  1.  The  plasmaforming  gas  mass  flow  rate  is  5  g/s,  the  swirling  angle  is  45°,  the  energy  input  is 
«15  kW,  viscosity  and  thermal  conductivity  effects  have  been  disregarded. 

B  2.  The  plasmaforming  gas  mass  flow  rate  is  2  g/s,  the  swirling  angle  is  45°,  the  energy  input  is 
«15  kW,  viscosity  and  thermal  conductivity  effects  have  been  disregarded. 

B3.  The  plasmaforming  gas  mass  flow  rate  is  2  g/s,  the  swirling  angle  is  45°,  the  energy  input  is 
«15  kW,  viscosity  and  thermal  conductivity  effects  have  been  accounted. 

B4.  The  plasmaforming  gas  mass  flow  rate  is  2  g/s,  the  swirling  angle  is  0°  (there  is  no  swirling), 
the  energy  input  is  «15  kW,  viscosity  and  thermal  conductivity  effects  have  been  accounted. 

B1  variant  In  Fig.3.Bla-c  one  can  see  distributions  of  the  pressure,  temperature  and  the 
velocity  absolute  value  (the  module  of  the  velocity)  in  the  conical  plasma  generator  at  the 
inputted  power  W~15  kW  (injection  to  the  counterflow).  The  plasmaforming  gas  mass  flow  rate 
is  5  g/s,  the  swirling  angle  is  45°,  viscosity  and  thermal  conductivity  effects  have  been 
disregarded. 
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Fig.3.Bla.  Pressure  distribution  in  the  conical  plasma  generator  at  the  inputted  power  of  W~15  kW  (injection  to 
the  counterflow).  The  plasma  forming  gas  mass  flow  rate  is  5  g/s,  the  swirling  angle  is  45°,  viscosity  and  thermal 
conductivity  effects  are  disregarded. 
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Fig.  3.Blb.  Temperature  distribution  in  the  conical  plasma  generator  at  the  inputted  power  of  W~15  kW 
(injection  to  the  counterflow).  The  plasma  forming  gas  mass  flow  rate  is  5  g/s,  the  swirling  angle  is  45°,  viscosity 
and  thermal  conductivity  effects  are  disregarded. 
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Fig.S.Blc.  The  velocity  module  distribution  in  the  conical  plasma  generator  at  the  inputted  power  of  W~15  kW 
(injection  to  the  counterflow).  The  plasma  forming  gas  mass  flow  rate  is  5  g/s,  the  swirling  angle  is  45°,  viscosity 
and  thermal  conductivity  effects  are  disregarded. 


B2  variant.  In  Fig.3.B2  a-c  one  can  see  distributions  of  the  pressure,  temperature  and  the 
velocity  absolute  value  (the  module  of  the  velocity)  in  the  conical  plasma  generator  at  the 
inputted  power  W~15  kW  (injection  to  the  counterflow).  The  plasma  forming  gas  mass  flow  rate 
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is  2  g/s,  the  swirling  angle  is  45°,  viscosity  and  thermal  conductivity  effects  have  been 
disregarded.. 
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Fig.3.B2a.  Pressure  distribution  in  the  conical  plasma  generator  at  the  inputted  power  of  W~15  kW  (injection  to  the 
counterflow).  The  plasma  forming  gas  mass  flow  rate  is  2  g/s,  the  swirling  angle  is  45°,  viscosity  and  thermal 
conductivity  effects  are  disregarded. 
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Fig.3.B2b.  temperature  distribution  in  the  conical  plasma  generator  at  the  inputted  power  of  W~15  kW  (injection  to 
the  counterflow).  The  plasma  forming  gas  mass  flow  rate  is  2  g/s,  the  swirling  angle  is  45°,  viscosity  and  thermal 
conductivity  effects  are  disregarded. 
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Fig.3.B2c.  Distribution  of  the  velocity  module  in  the  conical  plasma  generator  at  the  inputted  power  of  W~15  kW 
(injection  to  the  counterflow).  The  plasma  forming  gas  mass  flow  rate  is  2  g/s,  the  swirling  angle  is  45°,  viscosity 
and  thermal  conductivity  effects  are  disregarded. 


B3  variant.  In  Fig.3.B3  a-c  one  can  see  distributions  of  the  pressure,  temperature  and  the 
velocity  absolute  value  (the  module  of  the  velocity)  in  the  conical  plasma  generator  at  the 
inputted  power  W~15  kW  (injection  to  the  counterflow).  The  plasma  forming  gas  mass  flow  rate 
is  2  g/s,  the  swirling  angle  is  45°,  viscosity  and  thermal  conductivity  effects  have  been 
accounted. 
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Fig.3.B3a.  Pressure  distribution  in  the  conical  plasma  generator  at  the  inputted  power  of  W~15  kW  (injection  to  the 
counterflow).  The  plasma  forming  gas  mass  flow  rate  is  2  g/s,  the  swirling  angle  is  45°,  viscosity  and  thermal 
conductivity  effects  are  accounted. 
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Fig.3.B3b.  Temperature  distribution  in  the  conical  plasma  generator  at  the  inputted  power  of  W-15  kW  (injection  to 
the  counterflow).  The  plasma  forming  gas  mass  flow  rate  is  2  g/s,  the  swirling  angle  is  45°,  viscosity  and  thermal 
conductivity  effects  are  accounted. 
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Fig.3.B3c.  Distribution  of  the  velocity  module  in  the  conical  plasma  generator  at  the  inputted  power  of  W~15  kW 
(injection  to  the  counterflow).  The  plasma  forming  gas  mass  flow  rate  is  2  g/s,  the  swirling  angle  is  45°,  viscosity 
and  thermal  conductivity  effects  are  accounted. 


B4  variant.  In  Fig.3.B4  a-c  one  can  see  distributions  of  the  pressure,  temperature  and  the 
velocity  absolute  value  (the  module  of  the  velocity)  in  the  conical  plasma  generator  at  the 
inputted  power  W-15  kW  (injection  to  the  counterflow).  The  plasma  forming  gas  mass  flow  rate 
is  2  g/s,  the  swirling  angle  is  0°  (no  swirling),  viscosity  and  thermal  conductivity  effects  have 
been  accounted. 
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Fig.3.B4a.  Pressure  distribution  in  the  conical  plasma  generator  at  the  inputted  power  of  W-15  kW  (injection  to  the 
counterflow).  The  plasma  forming  gas  mass  flow  rate  is  2  g/s,  the  swirling  angle  is  0°,  viscosity  and  thermal 
conductivity  effects  are  accounted. 
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Fig.3.B4b.  Temperature  distribution  in  the  conical  plasma  generator  at  the  inputted  power  of  W~15  kW  (injection  to 
the  counterflow).  The  plasma  forming  gas  mass  flow  rate  is  2  g/s,  the  swirling  angle  is  0°,  viscosity  and  thermal 
conductivity  effects  are  accounted. 
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Fig.3.B3c.  Distribution  of  the  velocity  module  in  the  conical  plasma  generator  at  the  inputted  power  of  W~15  kW 
(injection  to  the  counterflow).  The  plasma  forming  gas  mass  flow  rate  is  2  g/s,  the  swirling  angle  is  0°,  viscosity 
and  thermal  conductivity  effects  are  accounted. 


The  comparison  of  B1  and  B2  computation  results  shows  that  the  increasing  of  the  mass 
flow  rate  of  the  plasma  forming  gas  at  other  same  parameters  leads  to  the  enlarging  and  mainly 
to  the  elongation  of  the  hot  area  outside  the  plasma  generator. 

The  comparison  of  B2  and  B3  computation  results  shows  that  accounting  of  thermal 
conductivity  and  viscosity  substantially  influences  distributions  of  temperature  and  the  absolute 
velocity  value  and  to  increasing  (widening  in  the  transversal  direction  by~1.5-2  times)  of  the 
region  with  high  values  of  parameters.  This  substantially  changes  the  temperature  distribution  in 
the  counterflow  outside  the  plasma  generator. 
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The  comparison  of  B3  and  B4  computation  results  shows  that  variation  of  the  swirling 
angle  weakly  influences  the  parameter  distribution  inside  and  outside  the  plasma  generator.  This 
says  that  the  swirling  is  necessary  only  for  the  stabilizing  of  the  jet  inside  the  plasma  generator 
and  to  prevent  melting  of  its  surface  due  to  the  development  of  different  instabilities  developed 
at  plasma  local  thermal  non-uniformity  what  is  correlated  with  the  results  of  [3.8].  This  result 
allows  us  to  use  parameter  distribution  computation  results  obtained  with  accounting  of  swirling 
and  without  of  it  at  the  investigation  of  the  plasma  jet  interaction  with  the  cross  supersonic  flow. 
The  last  conclusion  requires  additional  experimental  validation. 

3.6.  Interaction  of  flat  plasma  jet  containing  of  a  row  of  plasma  generators  with  conical 
nozzle. 

In  frames  of  works  on  the  project  we  carried  out  calculations  of  flat  plasma  jet  containing 
of  a  row  of  plasma  generators  with  conical  nozzle  and  creating  a  slot  (see  Chapter  2) 

Calculations  of  the  plasma  jet  from  the  conical  plasma  generator  (with  the  diverging 
channel)  and  the  supersonic  crossflow  interaction  at  angles  90°  and  135°  have  been  made  for 
determining  of  optimal  conditions  for  the  fuel  activation.  These  computations  are  of  the  principle 
interest  from  the  point  of  view  of  the  conical  plasma  generator  application  for  problems  of  the 
fuel  ignition  and  combustion.  In  this  case  we  have  the  comparison  opportunity  with  our 
computation  results  for  the  plasma  generator  with  the  cylindrical  channel  with  a  cavity  for 
plasma  mixing  ,  see  Chapter  2. 

The  flow  parameters  were  the  same  as  in  Chapter  2  and  section  3.4:  the  crossflow 
velocity  was  M=2,  the  pressure  was  0.2  atm,  the  temperature  was  250  K. 

As  earlier  the  mass  flow  rate  was  given  at  the  inflow  to  the  plasma  generator,  the  pressure 
is  extrapolated  inside  the  plasma  generator  to  the  inflow  boundary,  density  and  the  temperature 
are  determined  from  the  condition  of  the  isentropic  gas  transition  from  the  normal  state  to  the 
state  with  the  determined  pressure,  the  velocity  of  the  inflow  is  determined  by  the  given  mass 
flow  rate. 

In  the  present  calculations  we  use  the  data  of  B4  (section  3.4)  computations  on  the  jet 
radial  parameters  at  the  outlet  from  the  nozzle  when  the  gas  mass  flow  rate  was  2  g/s,  the 
swirling  angle  was  0°,  and  the  viscosity  and  thermal  conductivity  effects  were  accounted. 

However,  the  viscosity  and  thermal  conductivity  effects  at  computations  of  the  jet 
interaction  with  the  crossflow  outside  the  plasma  generator  proves  to  be  small  at  typical  times  of 
the  hydrodynamic  interaction  at  typical  gas  temperatures  outside  the  plasma  generator. 

In  Fig.3.5.1.  a-d  one  can  see  the  pressure,  temperature,  Mach  number  and  the  velocity 
absolute  value  (the  module  of  the  velocity)  distributions  at  the  plasma  jet  interaction  with  the 
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crossflow  at  90°  angle,  at  the  inputted  power  W~15  kW.  The  mass  flow  rate  of  the  plasma 
forming  gas  is  2  g/s,  the  swirling  angle  is  0°. 
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Fig. 3. 5.1. a  Pressure  distribution  at  the  plasma  jet  interaction  with  the  crossflow  at  90°  angle,  at  the  inputted  power 
ofW^15kW. 
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Fig.3.5.1.b  Temperature  distribution  at  the  plasma  jet  interaetion  with  the  crossflow  at  90°  angle,  at  the  inputted 
powerofW~15  kW. 
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Fig.3.5.1.c  Mach  number  distribution  at  the  plasma  jet  interaction  with  the  crossflow  at  90°  angle,  at  the  inputted 
power  of  W~15  kW. 
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Fig.3.5.1.d  The  module  of  the  velocity  distribution  at  the  plasma  jet  interaction  with  the  crossflow  at  90°  angle,  at 
the  inputted  power  of  W~15  kW. 


In  Fig.3.5.2.  a-d  one  can  see  the  pressure,  temperature,  Mach  number  and  the  velocity  absolute 
value  (the  module  of  the  velocity)  distributions  at  the  plasma  jet  interaction  with  the  crossflow  at 
135°  angle,  at  the  inputted  power  W~15  kW.  The  mass  flow  rate  of  the  plasma  forming  gas  is  2 
g/s,  the  swirling  angle  in  the  plasma  generator  is  0°. 

The  appearance  of  the  interaction  between  the  jet  and  the  flow  at  the  angle  90°  is 
practically  the  same  as  in  the  case  of  the  plasma  generator  with  the  cylindrical  channel  (with  the 
cavity)  analyzed  by  us  in  Chapter  2.  The  external  supersonic  crossflow  drifts  and  turns  the  jet 
down  the  flow.  In  the  case  of  Chapter  2  the  character  of  the  interaction  at  90°  and  135°  with 
crossflow  was  analogous. 

However,  the  comparison  of  interaction  of  the  jet  and  the  crossflow  in  case  of  the  conical 
plasma  generator  at  the  angles  90°  and  135°  reveals  the  substantial  differences. 
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Fig. 3. 5.2. a  Pressure  distribution  at  the  plasma  jet  interaction  with  the  crossflow  at  135°  angle,  at  the  inputted  power 
of  W~15  kW. 
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Fig.3.5.2.b  Temperature  distribution  at  the  plasma  jet  interaction  with  the  crossflow  at  135°  angle,  at  the  inputted 
power  ofW^15kW. 
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Fig. 3. 5.2. c  Mach  number  distribution  at  the  plasma  jet  interaction  with  the  crossflow  at  135°  angle,  at  the  inputted 
power  of  W~15  kW. 


Fig. 3. 5.2. d  The  module  of  the  velocity  distribution  at  the  plasma  jet  interaction  with  the  crossflow  at  135°  angle,  at 
the  inputted  power  of  W~15  kW. 
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Fig.  3. 5.2. a  demonstrates  that  deep  penetration  of  the  jet  towards  the  flow  takes  place.  Fig. 
3.5.2.b  demonstrates  that  the  gas  temperature  distribution  at  the  plasma  jet  interaction  with  the 
crossflow  at  135°  angle  is  extremely  non-uniform.  The  hot  region  is  brought  forth  towards  the 
flow.  So  the  region  of  the  plasma-flow  increases  in  comparison  with  the  case  of  90°.  Its 
temperature  is  3000  K.  It  is  quite  enough  for  the  ignition  of  the  hydrocarbon  fuels.  This  effect 
can  be  extremely  useful  in  problems  of  the  fuel  ignition  and  combustion.  Fig.  3.5.2.d 
demonstrates  that  the  incident  gas  flows  around  the  plasma  area.  At  that  the  interaction  time 
increases  in  comparison  with  the  case  of  90°,  and  it  is  in  its  turn  is  useful  at  initiating  of  chemical 
reactions. 

Undertook  computations  have  shown  the  prospects  of  the  conical  plasma  generator  with 
the  diverging  plasma  channel  application  for  the  problems  of  the  fuel  activation  and  ignition  in 
the  supersonic  flow. 

3.7.  Conclusions  to  Chapter  3 

Numerical  modeling  of  interaction  of  plasma  generator  with  divergent  nozzle  with 
motionless  gas  and  the  counterflow  was  made.  Comparison  of  computation  results  with  known 
experiments  showed  satisfactory  agreement.  It  allows  to  use  developed  approach  for  umodeling 
of  interaction  with  cross  flows  of  gases  and  flammable  mixture. 

The  comparison  of  computation  results  shows  that: 
increasing  of  the  mass  flow  rate  of  the  plasma  forming  gas  at  other  same  parameters  leads  to  the 
enlarging  and  mainly  to  the  elongation  of  the  hot  area  outside  the  plasma  generator; 
accounting  of  thermal  conductivity  and  viscosity  substantially  influences  distributions  of 
temperature  and  the  absolute  velocity  value  and  to  increasing  (widening  in  the  transversal 
direction  by~1.5-2  times)  of  the  region  with  high  values  of  parameters,  this  substantially  changes 
the  temperature  distribution  in  the  counterflow  outside  the  plasma  generator; 
variation  of  the  swirling  angle  weakly  influences  the  parameter  distribution  inside  and  outside 
the  plasma  generator.  This  says  that  the  swirling  is  necessary  only  for  the  stabilizing  of  the  jet 
inside  the  plasma  generator  and  to  prevent  melting  of  its  surface  due  to  the  development  of 
different  instabilities  developed  at  plasma  local  thermal  non-uniformity  what  is  correlated  with 
the  results  of  [3.8].  This  result  allows  us  to  use  parameter  distribution  computation  results 
obtained  with  accounting  of  swirling  and  without  of  it  at  the  investigation  of  the  plasma  jet 
interaction  with  the  cross  supersonic  flow.  The  last  conclusion  requires  additional  experimental 
validation. 

We  carried  out  calculations  of  flat  plasma  jet  containing  of  a  row  of  plasma  generators 
with  conical  nozzle  and  creating  a  slot  (see  Chapter  2)  interaction  with  the  cross  supersonic  flow 
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at  angles  90°  and  135°  for  possible  experimental  eonditions  for  determination  of  optimal 
eonditions  for  fuel  aetivation. 

Comparison  with  results  for  plasma  generator  with  eylindrieal  ehannel  shows  longer 
range  and  larger  area  of  plasma  jet  interaetion  with  the  eross  and  eounter  flow  in  the  ease  of 
plasma  generator  with  eonieal  nozzle.  Most  interesting  results  were  obtained  at  plasma  jet 
interaetion  with  the  eross  flow  at  the  angle  of  135°  .  The  gas  temperature  distribution  at  the 
plasma  jet  interaetion  with  the  erossflow  is  extremely  non-uniform  in  this  ease.  The  hot  region  is 
brought  forth  towards  the  flow.  So  the  region  of  the  plasma-flow  inereases  in  eomparison  with 
the  ease  of  90°.  Its  temperature  is  3000  K.  It  is  quite  enough  for  the  ignition  of  the  hydroearbon 
fuels.  This  effeet  ean  be  extremely  useful  in  problems  of  the  fuel  ignition  and  eombustion. 
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Chapter  4.  Interaction  of  nitrogen  plasma  jet  with  propane  air  mixture 


4.1.  Calculations  of  plasma  parameters  for  possible  experimental  conditions  when  nitrogen 
plasma  jet  is  injected  from  plasma  generator  to  motionless  propane  air  mixture. 

This  part  of  work  is  necessary  for  possible  experimental  modeling  and  choice  of  optimal 
hydrodynamic  conditions.  In  this  investigation  part  we  posed  the  problem  of  ignition  of 
stoichiometric  propane-air  mixture  by  the  nitrogen  plasma  jet  with  the  parameters:  plasma 
forming  mass  flow  rate  2  g/s,  temperature  of  the  gas  1000  K,  power  putted  to  plasma  is  5  kW  at 
the  gas  pressure  0.2  atm. 

The  solution  of  this  problem  was  necessary  before  modeling  of  plasma  jet- incident  flow  and  of 
plasma  jet-crossflow  interactions  for  choice  of  initial  conditions  of  plasma  jet  injection  to 
propane-air  mixture  at  some  angle  to  a  flow. 

Modeling  of  chemical  kinetics  of  the  stoichimetric  propane-air  mix  combustion. 

For  analysis  of  ignition  and  combustion  properties  of  the  hydrocarbons  in  air  it  is  necessary  to 
make  choice  of  main  chemical  reaction  leading  to  them.  A  number  of  works  on  the  propane 
combustion  is  rather  large  and  it  includes  different  approaches  to  the  combustion  modeling  [4.1  - 
7]  and  different  sets  of  chemical  reactions  taking  part  in  these  processes.  However  in  solution  of 
problems  connected  with  modeling  of  complex  gasdynamic  flows  together  with  combustion 
processes  it  is  necessary  reduce  as  possible  the  number  of  chemical  reagents  leaving  the  main 
features  of  the  processes.  In  this  direction  we  know  a  number  of  references  [4.8-1 1]  devoted  to 
development  and  application  of  the  simplified  schemes  for  description  of  complex  gasdynamic 
processes  with  combustion.  Using  results  of  works  [4.9-10]  in  our  investigations  of  gas 
discharges  in  the  propane-air  mix  [4.12-13]  we  obtained  times  of  the  mix  ignition  consistent  with 
our  experiments.  This  determined  the  choice  of  the  reference  [4.9]  as  starting  one  for  ignition 
and  combustion  of  stoichiometric  propane-air  mix  by  the  plasma  jet.  We  do  not  know  now  any 
work  devoted  to  ignition  of  flammable  mixes  by  plasma  jets.  So  there  was  no  opportunity  to 
compare  our  results. 

The  following  one  step  system  of  chemical  reactions  describes  well  variation  of  pressure  and 
component  composition  at  the  combustion  of  the  propane-air  stoichiometric  propane  air  mix 
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nl  and  n2  correspond  to  weight  eoneentrations  of  propane  and  oxygen,  u  is  the  volume  energy 
density  released  at  combustion.  Coefficients  in  the  equations  above  have  the  following  form  : 
cl  =  ^4  •  exp(—E  /T ) 

kl  =  (ml •  (m2 •  cl  •  (nl  f  ■  (nl / 

k2  =  5.Q-  (ml •  (ml •  cl  •  (nl  f  ■  (nl /  (4. 1 .4) 

k3  =  (ml •  (ml )-^  •  cl  •  (nl ^  •  (nl /  •  AH 
A=10'^;  E=15  100  K[4.9]; 

T  is  in  K,  ml,  m2  are  the  molar  masses  of  the  propane  and  oxygen  respectively  (at  normal 
eonditions  they  are  ml=44,  m2=32),  a=0.1;  P=1.65  [4.9];  the  dimensions:  [ni]=[g/cm  ],  and  [u] 
=[J/em^],  so  the  right  hand  sides  of  equations  (4.1.1)-(4.1.2)  are  expressed  in  g/(cm^  s),  and  the 
right  hand  sides  of  the  equation  (4. 1 .3)  is  expressed  in  J/(em^  s).  The  propane  -air  mix  has  the 
following  eomponent  composition  N2:  O2:  C3H8=0.7465:  0.2105:  0.043,  it  was  determined  by  the 
standard  equation  [4.14-15] 

C3//3  +  5  O2  +  5  •  3.762  3CO^  +  4H^O  +  5  •  3.762 

The  formation  of  substances  enthalpy  AH  was  ealculated  by  us  on  a  basis  of  the  standard 
equation  [4.14-15]: 

=  Ac./,.  (T)  +  5h„^  (T)  -  3hco,  (T)  -  4h^,o(T) 

with  a  help  of  data  represented  in  [4.16],  and  it  proved  to  be  AH  =  (2.0810.04)- 10^  [J/mol]  in  the 
range  200-4000  K,  whieh  corresponds  to  the  temperature  range  neeessary  at  ignition  and 
combustion  of  propane  -air  mixes 

4.1.2  Modeling  of  the  plasma  gas  dynamics  of  the  propane  -air  mix  combustion. 

At  theoretical  and  numerical  modeling  we  used  the  model  based  on  the  following 
equation  system  of  radiative  plasma  eleetrodynamies  represented  in  Chapter  1 . 

In  the  present  case  stoichiometric  propane-air  mix  parameters,  enthalpy  of  substances 
formation  and  ehemieal  reaetions  of  propane  and  oxygen  were  given  outside  the  plasma 
generator  and  plasma  jet,  and  inside  them  the  problem  was  solved  as  it  was  in  Chapters  2-3. 

For  eomputations  we  used  Lagrange  method  on  a  completely  eonservative  implieit 
differenee  scheme  with  an  adaptive  triangular  unstructured  grid.  This  was  effeetive  in  modeling 
the  eoupling  of  the  internal  and  external  flows  inside  and  outside  the  counterflow  jet  and  its 
nozzle  ,  see  Chapter  2-3.  The  total  number  of  knots  was  15000-150000.  The  grid  eoneentrated  in 
the  region  of  the  combustion  by  the  adaptive  method. 
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Chapter  4.1.3  Results  of  plasma  jet  from  the  plasma  generator  injection  to 
motionless  stoichiometric  mixture. 

Obtained  data  at  the  gas  pressure  0.2  atm  ,  power  putted  in  the  gas  5  kW,  initial  mixture 
temperature  T=1000  K  and  plasma  forming  gas  mass  flow  rate  2  g/s  are  represented  in  Fig. 
4.1.3. 1-4.1.3.5. 
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Fig.  4.1.3. 1  Pressure  distribution  at  plasma  jet  injection  to  stoichiometric  propane  air 


mixture  at  time  moment  0.5  •10"'^  s.  Pressure  0.2  atm  ,  Power  5  kW,  initial  temperature  T=1000  K. 


Fig.  4.1.3.2  Mach  number  (mach=  |v|/|a  (local)  I )  distribution  at  plasma  jet  injection  to 
stoichiometric  propane  air  mixture  at  time  moment  0.5  lO  "^  s.  Pressure  0.2  atm  ,  Power  5  kW,  initial 


temperature  T=1000  K. 


54 


E 


o 


tiine=o.50E+(>!>  Total  density,  10  kg/m® 


1  58E-D1 
1 .52E-01 
1  4SE-0t 
1  39E-01 
1  33E-01 
1 .27E-01 
1 .20E-01 
1,14E  01 
1  .OSE-01 
1.01E-0t 
9.50E-02 
S86E-02 
8  23E-a2 
7  60E-02 
6  97E  02 
6  33E-02 
5  70E-02 
5.07E-02 
4  44E-02 
3.80E-02 
3.17E^02 
2,54E-02 
1.9  IE-02 
1 .27E-02 
6.41  E-03 


Fig.  4.1.3.3  Density  distribution  at  plasma  jet  injection  to  stoichiometric  propane  air  mixture  at 
time  moment  0.5  s.  Pressure  0.2  atm  ,  Power  5  kW,  initial  temperature  T=1000  K. 
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Fig.  4.1.3.4  Temperature  distribution  at  plasma  jet  injection  to  stoichiometric  propane  air  mixture 


at  time  moment  0.5  -lO'"^  s.  Pressure  0.2  atm  ,  Power  5  kW,  initial  temperature  T=1000  K. 
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Fig.  4.1.3.5  Propane  density  distribution  at  plasma  jet  injection  to  stoichiometric  propane  air 
mixture  at  time  moment  0.5  -lO'"^  s.  Pressure  0.2  atm  ,  Power  5  kW,  initial  temperature  T=1000  K. 
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The  given  power  of  5  kW  putted  in  the  gas  was  chosen  basing  on  results  of  Chapter  3,  as  the 
power  at  which  the  ignition  of  the  flammable  mixture  can  take  place. 

Analysis  of  results  represented  in  Fig.  4. 1.3. 1-4. 1.3. 5.  shows  that  the  presence  of  ignition 
source  over  the  jet  surface  leads  to  the  reconstruction  of  gasdynamic  flows.  In  particular  pressure 
and  mixture  density  around  the  jet  are  reconstructed.  The  ignition  takes  place  only  in  thin  layer 
near  the  side  jet  surface.  At  this  power  and  at  the  considered  computation  times  there  is  not  the 
complete  combustion  of  the  fuel. 

However,  our  calculations  show  that  the  jet  at  the  given  parameters  can  be  applied  for  the 
ignition  of  the  flammable  mixture. 

4.2.1  Calculation  of  flow  parameters  when  the  injection  of  cylindrical  plasma  jet  takes 
place  to  the  incident  supersonic  flow  of  stoichiometric  propane-air  mixture. 

This  part  of  work  was  necessary  from  the  point  of  view  of  modeling  and  choice  of  optimal 
hydrodynamic  conditions. 

In  this  part  of  investigations  we  posed  the  problem  of  the  stoichiometric  mixture  propane- 
air  ignition  oncoming  with  Mach  number  M=2  on  the  plasma  jet  with  parameters:  plasma 
forming  gas  mass  flow  rate  2  g/s,  gas  temperature  of  the  oncoming  flow  1000  K  at  pressure  0.2 
atm.  Jet  parameters  corresponded  to  those  investigated  by  us  in  Chapters  2-3,  when  we 
considered  the  jet  interaction  with  the  incident  air.  The  plasma  generator  power  was  15  kW  since 
in  previous  calculations  we  clarified  that  the  plasma  penetration  of  the  incident  flow  takes  place 
at  the  power  putted  in  the  gas  of  >  13  kW. 

The  solution  of  this  problem  was  necessary  before  modeling  of  plasma  jet-incident  flow  and  of 
plasma  jet-crossflow  interactions  for  choice  of  initial  conditions  of  plasma  jet  injection  to 
propane-air  mixture  at  some  angle  to  a  flow. 

The  general  system  of  equations  added  with  the  balance  equations  of  oxygen  and  propane  in 
represented  above.  The  calculations  were  made  with  the  adaptive  grid  in  the  combustion  region, 
the  flow  at  T=1000  K  comes  on  the  get,  the  pressure  is  0.2  atm.  The  plasma  generator  power  is 
15  kW,  the  plasma  forming  gas  mass  flow  rate  is  2  g/s. 

In  Fig.  4.2.1 -4.2.4  one  can  see  results  of  calculations  for  main  gasdynamic  parameters  realizing 
at  the  plasma  jet  interaction  with  the  incident  flow  of  the  flammable  (stoichiometric  propane-air) 
mixture. 

The  calculations  show  that  the  incident  flow  stops  and  turns  back  the  plasma  jet.  The  shock  wave 
appears  at  that,  pressure  and  temperature  rise  in  the  area  of  the  interaction,  and  the  combustion 
wave  is  realized.  It  propagates  in  all  directions  from  the  jet.  At  its  reflection  from  the  limiting 
wall  of  the  volume  there  appears  a  typical  Mach  foot.  At  that  the  temperature  in  the  air  of  gas 
dynamic  interaction  rises  up  to  4000  K,  caused  by  the  propane  combustion. 
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In  Fig.4.2.5  one  can  see  results  of  propane  density  calculations  that  is  realized  at  the  plasma  jet 
interaction  with  the  incident  hot  flow.  One  can  see  that  there  is  a  rise  of  propane  density  in  the 
front,  it  leads  to  acceleration  of  the  combustion  process. 

Analysis  of  calculated  data  shows  that  at  the  interaction  of  the  jet  with  propane-air  mix  there 
takes  place  the  ignition  of  the  mix  and  the  detonation  realization  during  the  time  of  ~0.6  10'"^  s. 
This  reveals  the  availability  of  pulse  periodic  plasma  jets  for  realization  of  the  detonation 
ignition  of  the  flammable  mixtures. 
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Fig.4.2.  1.  Plasma  jet  injection  to  stoichiometric  propane-air  mixture:  pressure  distribution.  Initial  pressure  in 
the  incident  flow  is  0.2  atm,  initial  gas  temperature  is  1000  K,  plasma  generator  power  is  15  kW.  Computation  time 
is  1.1  lO'^s. 
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Fig.4.2.  2.  Plasma  jet  injection  to  stoichiometric  propane-air  mixture:  total  density  distribution.  Initial 
pressure  in  the  incident  flow  is  0.2  atm,  initial  gas  temperature  is  1000  K,  plasma  generator  power  is  15  kW. 
Computation  time  is  l.llO'^s. 
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Fig.4.2.3.  Plasma  jet  injection  to  stoichiometric  propane-air  mixture:  temperature  distribution.  Initial 
pressure  in  the  incident  flow  is  0.2  atm,  initial  gas  temperature  is  1000  K,  plasma  generator  power  is  15  kW. 
Computation  time  is  1.1  10"^  s. 
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Fig.4.2.  4.  Plasma  jet  injection  to  stoichiometric  propane-air  mixture:  Mach  number  distribution.  Initial 
pressure  in  the  incident  flow  is  0.2  atm,  initial  gas  temperature  is  1000  K,  plasma  generator  power  is  15  kW. 
Computation  time  is  1.1  10"^  s. 
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Fig.4.2.5.  Plasma  jet  injection  to  stoichiometric  propane-air  mixture:  temperature  distribution.  Initial  pressure 
in  the  incident  flow  is  0.2  atm,  initial  gas  temperature  is  1000  K,  plasma  generator  power  is  15  kW.  Computation 
time  is  1.110"^  s. 

4.3.Parameters  of  flow  at  the  injection  of  the  flat  nitrogen  plasma  jet  from  the  plasma 
generator  to  crossflow. 

In  Fig.  4.3.1.  we  remind  the  displaeement  of  the  jet  (a  row  of  plasma  generators)  and  the  flow  at 
whieh  we  eonsider  their  interaetion. 


Fig.  4.3.1.  Flat  plasma  jet  created  by  a  row  of  plasma  generators 


The  present  eomputations  on  the  adaptive  grid  have  the  prineiple  eharaeter  for 
elarifieation  of  the  interaetion  features  between  the  plasma  jet  and  the  eross  flow.  The  have  been 
earned  out  to  realize  high  resolution  and  high  aeeuraey  in  the  region  of  interaetion.  Previous 
ealeulations  in  Capter  2-3,  were  earried  out  with  grid,  whieh  was  rather  detailed  only  near  the 
boundary  of  the  plasma  blow  in  the  erossflow.  The  grid  was  rough  far  from  this  boundary  and 
did  not  allow  to  aeeurately  determine  the  jet  and  flow  parameters  in  the  region  of  their 
interaetion. 

Method  of  the  iet  boundary  ealeulation 

The  ealeulation  is  made  by  Free  Lagrange  method  on  the  triangular  mesh:  the  ealeulated 
points  are  moving  together  with  the  medium.  The  ehanging  of  the  mesh  strueture  and  mesh 
adaptation  are  made  at  eaeh  step.  The  adaptation  is  undertaken  so  that  there  were  suffieient 
ealeulated  points  in  the  jet  and  in  the  flow  exeited  by  it  for  good  reproduetion  of  the  flow 


strueture. 
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The  blowing  of  the  jet  in  the  numerical  method  is  realized  by  the  following  way.  The  gas 
velocity  is  given  in  the  boundary  points  (lying  on  the  boundary)  of  the  numerical  domain  from 
which  the  blowing  takes  place.  The  thermodynamic  parameters  of  the  jet  are  given  in  the  near 
boundary  cells.  These  points  are  moving  with  the  given  velocity  into  the  numerical  domain  at 
each  time  step.  In  the  beginning  of  the  next  step  these  points  return  to  the  initial  position  at  the 
boundary  from  which  the  blowing  takes  place.  Such  a  procedure  of  the  blowing  realizing  in  the 
Lagrange  coordinates  approximates  well  the  boundary  conditions  of  the  blowing.  The  distances 
between  the  boundary  and  near  boundary  points  increase  in  the  result  of  this  procedure.  When 
these  distances  exceed  the  typical  quantity  of  the  calculation  mesh  step  then  the  new  nodes  are 
added  between  the  boundary  and  near  boundary  nodes. 

In  order  to  back  trace  the  jet  boundary  a  new  computation  variable  -  an  indication  is 
introduced,  it  is  specified  in  all  calculated  nodes.  The  indication  is  equal  to  zero  in  all  the  points 
of  the  incident  flow.  The  indication  is  equal  to  unity  in  the  points  lying  on  the  boundary  from 
which  the  blowing  takes  place.  At  addition  of  a  new  node  to  the  mesh  edge  the  indication  in  a 
new  point  is  determined  by  the  indications  at  the  edge  ends  by  linear  interpolation.  So  the 
indication  is  larger  then  zero  at  the  addition  of  new  nodes  at  the  blowing  boundary.  With  time 
the  indication  of  adding  nodes  after  multiple  additions  of  new  nodes  will  be  close  to  unity  at  the 
boundary  of  blowing.  And  the  indication  will  be  more  than  zero  and  close  to  unity  in  the 
calculated  points  corresponding  to  the  blowing  jet.  In  figures  the  boundary  of  the  jet  is  depicted 
by  the  indication  level  line  corresponding  to  its  definite  value,  which  is  close  to  unity.  In  our  case 
this  value  is  equal  to  0.8.  Values  of  the  indication  function  lies  in  the  range  between  0  and  1,  our 
calculations  show  that  all  the  indication  level  lines  (between  0  and  1)  lie  in  narrow  region 
approximately  determining  the  boundary  between  the  blowing  jet  and  the  incident  flow. 

Results  on  the  flat  jet,  (interaction  with  a  supersonic  flow). 

Gas  flow  parameters:  gas  temperature  is  T=250  K,  pressure  is  0.2  atm.,  Mach  number  is 

2. 

Plasma  jet  parameters:  nitrogen  jet  comes  out  of  the  plasma  generator,  electric  source 
power  P=  10  kW,  gas  flow  rate  2  g/s,  plasma  is  directed  at  the  angle  of  135°  to  the  cross  flow. 

In  Fig.  4.3.2-  4. 3. Tone  can  see  results  of  calculations  with  indication  of  the  boundary 
between  the  jet  and  the  flow  for  the  time  moment  t  =  3  - 10'^  cex. 

Calculations  show  that  the  jet  creates  a  shock  wave  in  the  gas  of  the  crossflow.  In 
particular  the  temperature  in  the  gas  of  the  cross  flow  exceeds  1500°  K  in  the  head  part.  It  is 
probable  that  the  heating  in  the  shear  boundary  layer  is  sufficient  for  the  fuel  ignition. 
Examination  of  this  idea  requires  additional  calculations  with  modeling  of  processes  in  the  shear 
boundary  layer. 
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It  is  necessary  to  emphasize  that  the  plasma  generator’s  jet  is  strongly  non-imifomi.  Its 
external  layers  interacting  with  air  crossflow  are  relatively  cold  (  their  temperatuie  T~  500°  K) 
and  do  not  contain  the  charged  plasma.  Internal  jet  layers  consist  of  such  a  plasma,  they  do  not 
interact  with  ah  flow.  The  jet  gas  expands  and  cools  to  temperatuie  of  T~  1000°  -  2000°  K  at 
moving  away  fiom  the  bomidary  of  the  blowing  in.  The  jet  in  this  case  consists  of  the  gas  and 
does  not  contam  the  charged  plasma  (see.  Fig.  4.3.7).  So  such  a  plasma  stays  in  the  jet  near  the 
bomidary  of  the  blowing  in  and  it  does  not  interact  with  the  air  cross  flow.  At  the  same  time 
there  is  practically  full  dissociation  of  nitrogen  molecules  m  the  plasma  near  the  region  of  the  jet 
injection  to  a  flow  (see  Fig.  4.3.8  for  atoms  Nn).  This  could  lead  to  additional  heating  of  this 
region,  due  to  diffiision  of  atoms  from  the  plasma  and  the  following  recombination  of  atoms. 

hi  figmes  one  can  see  that  the  jet  is  unstable,  its  heavy  periphery  parts  strongly  slow 
down  by  the  flow  and  can  mix  with  air.  However  the  temperature  in  the  mixing  zone  does  not 
exceed  T~1000-2000°  K.  This  temperatme  range  also  satisfies  conditions  of  hydrocarbon-air 
mixtures  ignition. 


Fig.  4.3.2  Distributiou  of  pressure  at  plasma  jet  iuteractiom  Ttith  the  supersonic  cross  flow  at  the  angle  of  135°. 
Cross  flow  paianieters:  T=250  K,  pressure  0.2  atui.,  Mach  number  M=  2.  Nitrogen  plasma  jet  parameters  plasma 
generator  power  is  W~10  kW,  plasma  forming  gas  mass  f  low  rate  is  2  g/s. 
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Fig.  4.3.3  Distribution  of  densitj^  at  plasma  Jet  iuteraction  vdth  the  supersonic  cross  flow  at  the  angle  of  135®. 
Cross  flow  paiameters:  T=250  K,  pressure  0.2  atm.,  Mach  number  M=  2.  Nitrogen  plasma  jet  parameters  plasma 
generator  power  is  W--10  kW,  plasma  forming  gas  mass  f  low  rate  is  2  g/s. 
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Fig.  4.3.4  Mach  number  distribution  at  plasma  jet  interaction  with  the  supersonic  cross  flow  at  the  angle  of 
135°.  Cross  flow  parameters:  T=250  K,  pressure  0.2  atm.,  Mach  number  M=  2.  Nitrogen  plasma  jet  parameters 
plasma  generator  power  is  W~10  kW,  plasma  forming  gas  mass  flow  rate  is  2  g/s. 


Fig.  4.3..5  Temperature  distribution  at  plasma  jet  interaction  with  the  supersonic  cross  flow  at  the  angle  of 
135°.  Cross  flow  parameters:  T=250  K,  pressure  0.2  atm.,  Mach  number  M=  2.  Nitrogen  plasma  jet  parameters 
plasma  generator  power  is  W~10  kW,  plasma  forming  gas  mass  flow  rate  is  2  g/s. 
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Fig.  4.3.6  Absolute  value  of  velocity  distribution  at  plasma  jet  interaction  with  the  supersonic  cross  flow  at  the 
angle  of  135°.  Cross  flow  parameters:  T=250  K,  pressure  0.2  atm.,  Mach  number  M=  2.  Nitrogen  plasma  jet 
parameters  plasma  generator  power  is  W~10  kW,  plasma  forming  gas  mass  flow  rate  is  2  g/s. 
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Fig.  4.3.7  Electron  concentration  distribution  at  plasma  jet  interaction  with  the  supersonic  cross  flow  at  the 
angle  of  135°.  Cross  flow  parameters:  T=250  K,  pressure  0.2  atm.,  Mach  number  M=  2.  Nitrogen  plasma  jet 
parameters  plasma  generator  power  is  W-10  kW,  plasma  forming  gas  mass  flow  rate  is  2  g/s. 
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Fig.  4.3..8  Atoms  concentration  at  plasma  jet  interaction  with  the  supersonic  cross  flow  at  the  angle  of  135°. 

Cross  flow  parameters:  T=250  K,  pressure  0.2  atm.,  Mach  number  M=  2.  Nitrogen  plasma  jet  parameters  plasma 
generator  power  is  W-10  kW,  plasma  forming  gas  mass  flow  rate  is  2  g/s. 


4.4.  Modeling  of  propane-air  mixture  combustion  at  plasma  jent  and  the  cross  flow 
interaction  at  the  angle  of  150° 


As  in  our  simulations  in  Chapter  2-3  the  gasdynamie  equations  are  eoupled  with  the  2-D 
Maxwell  equations  and  the  plasma  dynamies  equations,  aeeounting  for  ohmie  heating  and 
radiative  plasma  eooling.  Speeifie  parameters  are  nitrogen  gas  mass  flow  in  the  generator  mf  =  2 
g/s,  power  input  P=15,  30  kW.  Appropriate  thermodynamie,  optieal  and  transport  properties  of 
nitrogen  and  air  plasmas  were  used  in  the  eomputations.  Propane-air  ehemistry  was  the  same  as 
in  seetion  4.2.  A  Lagrange  method  based  on  a  eompletely  eonservative  implieit  differenee 
seheme  with  an  adaptive  triangular  unstruetured  mesh  deseribed  above  was  applied.  This  was 
effeetive  in  modeling  the  eoupling  of  the  internal  and  external  flows  inside  and  outside  the 
eounterflow  jet  and  its  nozzle.  The  mesh  was  eoneentrated  in  the  region  of  the  eombustion  by 
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the  adaptive  method.  Computations  of  plasma  jet  mjection  to  the  crossflow  of  propane-ah 
flammable  mix  have  been  made  at  the  angle  150  Plasma  forming  gas  (nitrogen)  mass  flow  rate 
was  2  g/s,  plasma  generator  power  was  15  kW  and  30  kW  at  the  temperatiue  of  the  incident 
propane-air  stoichiometric  mix  1000  K. 

In  Fig.  4.4. 1-6.  represented  results  of  computations  of  gasdynamic  and  fuel  characteristics 
distributions  at  the  plasma  jet  interaction  with  the  supersonic  propane-ah  stoichiometric  mix 
flow  at  the  angle  150°  for  the  moment  18  ps.  The  flow  parameters:  T=1000  K,  Pressme  0.2  atm., 
the  Mach  number  is  M=  2.  The  jet  parameters:  the  nitrogen  jet  from  the  plasma  generator.  The 
soiu'ce  power  W=15  kW,  the  mass  flow  rate  2  g/s  (up  down  in  the  figmes). 
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Fig.4.4.1.  Pressure  distribution  at  flat  plasma  jet  interactiou  \i1tli  the  supersonic  propane-air  crossflow  (flow 
direction  is  up-down)  under  the  angle  of  150®  at  the  time  moment  18  ps.  Flow  parameters:  T=1 000  K, 
pressiu  e  0.2  atm.,  M=  2.  Nitiogeo  plasma  jet  parameters:  putted  inpow’^ei  is  W=15  kW,  plasma  foiming  gas  mass 
flow  rate  is  2  g/s. 
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Fig.4.4.2.  Mach  number  distribution  at  flat  plasma  jet  interaction  with  the  supersonic  propane-air  crossflow 
(flow  direction  is  up-down)  under  the  angle  of  150°  at  the  time  moment  18  ps.  Flow  parameters:  T=1000  K, 
pressure  0.2  atm.,  M=  2.  Nitrogen  plasma  jet  parameters:  putted  in  power  is  W=15  kW,  plasma  forming  gas  mass 
flow  rate  is  2  g/s. 
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Fig.4.4.3.  Total  density  distribution  at  flat  plasma  jet  interaction  with  the  supersonic  propane-air  crossflow 
(flow  direction  is  up-down)  under  the  angle  of  150°  at  the  time  moment  18  ps.  Flow  parameters:  T=1000  K, 
pressure  0.2  atm.,  M=  2.  Nitrogen  plasma  jet  parameters:  putted  in  power  is  W=15  kW,  plasma  forming  gas  mass 
flow  rate  is  2  g/s. 
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Fig.4.4.4.  Temperature  distribution  at  flat  plasma  jet  interaction  with  the  supersonic  propane-air  crossflow 
(flow  direction  is  up-down)  under  the  angle  of  150°  at  the  time  moment  18  ps.  Flow  parameters:  T=1000  K, 
pressure  0.2  atm.,  M=  2.  Nitrogen  plasma  jet  parameters:  putted  in  power  is  W=15  kW,  plasma  forming  gas  mass 
flow  rate  is  2  g/s. 
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Fig.4.4.5.  Oxygen  density  distribution  at  flat  plasma  jet  interaction  with  the  supersonic  propane-air 
crossflow  (flow  direction  is  up-down)  under  the  angle  of  150°  at  the  time  moment  18  ps.  Flow  parameters: 
T=1000  K,  pressure  0.2  atm.,  M=  2.  Nitrogen  plasma  jet  parameters:  putted  in  power  is  W=15  kW,  plasma  forming 
gas  mass  flow  rate  is  2  g/s. 
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Fig.4.4.6.  Fropaue  dens^it^-  distributiou  at  flat  plasma  Jet  interaction  \Yith  the  snpersonic  propane- air 
crossflow  (flow  direction  is  up-down)  under  the  angle  of  150®  at  the  time  moment  18  ps.  Flow’^  parameters: 
T=1000  K,  pressiue  0.2  atin.^  M=  2.  Nitiogen  plasma  jet  paiameters:  putted  in  power  is  W=15  kW,  plasma  forming 
gas  mass  flow  rate  is  2  g/s. 

Comparison  of  Figines  4.4.5-6  and  4.4.4  shows  quick  propagation  of  density,  pressure  and 
temperature  wave  in  case  of  the  fuel  mixtiue.  At  that  it  propagates  coimter  the  supersonic  flow, 
and  this  validates  the  occurrence  of  detonation.  At  the  same  time  considerable  increase  of  density 
and  pressiue  takes  place  on  the  wave  front,  it  leads  to  acceleration  of  propane  combustion.  In  the 
jet  area  appear  regions  with  low  density  it  indicates  appearance  of  vorticle  structures,  there  takes 
place  the  break  of  the  jet  also,  hi  the  result  there  appears  a  large  area  of  the  hydrodynamic 
interaction  of  the  plasma  with  the  crossflow.  Results  of  computations  represented  in  Fig.4.4.5 
and  4.4.6  show  fast  propane  burning  out  and  oxygen  decrease  in  the  region  of  high  temperabire, 
where  it  reaches  3000-4000  K.  Our  calculations  show  that  realization  of  combustion  wave  takes 
place  diuing  —  8-10  p,s. 

Computation  results  show  the  possibility  of  propane-air  mixtme  crossflow  ignition  at  its 
interaction  with  the  nitrogen  plasma  jet  from  the  plasma  generator  with  power  putted  into  the  gas 
of  P  =  15  kW  at  the  initial  temperature  of  1000  K. 

The  mixtuie  is  compressed  in  the  combustion  wave  front  (see  red  curve  on  the  front  in 
Fig  4.4.6)  and  it  accelerates  ignition. 

We  do  not  know  now  any  work  devoted  to  ignition  of  flammable  supersonic  mixes  by 
plasma  jets.  So  there  was  no  opportunity  to  compare  oui'  results.  Deteimination  of  the  minimum 
plasma  jet  power  at  which  the  plasma  from  the  jet  penetrates  the  cross  flow  and  still  ignites  the 
fuel  requites  additional  investigations.  However,  we  are  limited  with  power  P  >13  at  which 
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plasma  jet  is  locked  by  the  cross  flow.  Our  results  show  the  possibility  of  application  of  high 
power-low  energy  plasma  jets  with  the  pulse  duration  of  about  10  -100  ps  for  ignition  and 
combustion  applications.  They  can  be  used  in  pulse  repetitive  mode  which  decrease  the  energy 
losses  of  plasma  jets. 

Conclusions  to  Chapter  4 

1 .  The  model  of  ignition  and  combustion  of  the  stoichiometric  propane-air  mixture  with  a 
help  of  the  plasma  jet  from  the  plasma  generator  to  the  incident  flow  and  crossflow  of  the 
flammable  mix.  This  model  uses  one  step  effective  chemical  reaction  between  propane 
and  oxygen  molecules.  It  accounts  heating  of  the  medium  due  to  energy  released  by 
chemical  reactions  and  by  the  plasma  jet. 

2.  Lagrange  method  on  a  completely  conservative  implicit  difference  scheme  with  an 
adaptive  triangular  unstructured  grid  was  used  in  computations.  The  total  number  of 
knots  was  15000-150000.  The  mesh  was  concentrated  in  the  region  of  the  combustion  by 
the  adaptive  method. 

3.  This  method  is  applied  for  clarification  of  initial  and  boundary  conditions  at  investigation 
of  the  plasma  jet- flammable  mix  cross  flow  interaction  (at  angles  135,  150”). 

4.  Computations  of  plasma  jet  injection  to  the  incident  flow  of  propane-air  flammable  mix 
have  been  made.  Plasma  forming  gas  (nitrogen)  mass  flow  rate  was  2  g/s,  plasma 
generator  power  was  15  kW  and  30  kW  at  the  temperature  of  the  incident  propane-air 
stoichiometric  mix  250  K  and  1000  K. 

5.  Results  of  computations  demonstrated  the  ignition  possibility  of  propane-air  mixture  with 
a  help  of  plasma  generators  with  the  power  15  kW  at  the  incident  propane-air 
stoichiometric  mix  temperature  1000  K. 

6.  The  realization  of  fast  combustion  wave  processes  at  typical  times  of  ~10  microseconds 
is  observed  both  in  cylindrical  and  flat  cases. 
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Chapter  5.  Combined  discharge  plasma  influence  on  propane-air  mix  ignition 

5.1  Formulation.  Ignition  enhancement  of  fuels  with  a  help  of  plasma  is  under  detailed 
investigation  of  last  years.  Our  works  on  this  subject  devoted  to  an  application  of  plasma  jets 
have  shown  that  there  is  a  strong  gasdynamic  interaction  between  the  incident  flow  and  the  jet 
plasma  which  leads  to  fast  heating  of  the  incident  flammable  mix  with  creation  of  a  shock  in  the 
place  of  the  contact  and  the  following  heating  of  the  incident  mix  to  temperatures  necessary  for 
ignition.  However  plasma  chemical  reactions  do  not  play  important  role  due  to  slow  molecular 
transport  from  the  plasma  zone  to  the  fuel.  Analogous  situation  is  observed  in  some  works  on 
heat  sources  application  in  supersonic  combustion  [5.1].  In  this  case  it  is  reasonable  to  locate 
some  sort  of  plasma  source  to  insure  plasma  chemical  influence  on  the  fuel  with  production  of 
excited  species  and  active  radicals.  In  order  such  a  source  could  produce  active  radicals 
independently  of  gas  flow  parameters  the  time  of  its  interaction  has  to  be  smaller  than  typical 
gasdynamic  times. 

There  are  a  number  of  devices  which  can  produce  such  an  effect,  they  are  X-radiation 
plates,  non-selfmaintained  discharges,  and  ultra-violet  radiation  sources.  Application  of  X- 
radiation  requires  high  parameters  of  radiation  sources,  which  are  not  realized  now  and 
applications  of  special  protection  measures.  In  this  work  we  will  consider  an  application  of  the 
non-selfmaintained  discharge.  In  Fig.l  one  can  see  a  possible  test  channel  scheme  with  the 
application  of  the  non-selfmaintained  discharge  and  plasma  jet.  The  self  maintained  discharge 
device  is  located  upstream  the  step  and  the  plasma  jet.  The  self  maintained  discharge  represents  a 
combination  of  electron  beam  and  the  gas  discharge.  Where  electrons  and  ions  are  produced  by 
the  beam,  and  the  discharge  is  realized  at  voltage  lower  than  the  breakdown  one.  In  this  case 
produced  excited  components  will  be  moved  down  where  they  will  take  part  in  ignition  process. 

For  analysis  of  this  scheme  is  necessary:  a)  to  investigate  possible  parameters  of  the 
ignition  which  can  be  produced  by  the  flammable  mix  interaction  with  excited  chemically  active 
species  produced  by  the  discharge,  b)  to  formulate  necessary  parameters  for  non-selfmaintained 
discharge.  Namely  to  this  research  problem  this  work  is  devoted,  where  on  the  basis  of 
calculations  of  ignition  of  modified  stoichiometric  propane-air  mix  by  oxygen  atoms  generated 
by  non-selfmaintained  discharge  in  air  is  considered  the  question  of  applicability  of  this  scheme. 
5.2.  Approach.  For  solution  of  the  complete  problem  of  ignition  enhancement  of  propane-air 
stoichimetric  mixture  it  is  necessary  to  solve  the  problem  for  definite  channel  construction  with 
full-scale  solutions  of  gasdynamics  and  plasma  chemistry  equations  as  we  did,  for  example,  in 
[5.2].  Since  the  present  problem  has  the  research  character  we  consider  the  simplest  case  of 
addition  of  chemically  active  species  to  propane  air  stoichiometric  mix  in  motionless  conditions. 


70 


In  literature  [5.3-6]  there  is  an  information  only  about  reaction  of  the  atomic  oxygen  with 
components  of  propane  (CsHg)  decomposition,  and  indication  [5.6]  on  influence  of  vibrationally 
excited  state  of  the  oxygen  molecule  ,  02(v).  So  we  consider  the  ignition  of  propane-air 
stoichiometric  mix  at  addition  to  it  of  oxygen  atoms. 

The  model  of  CsHg  -air  mix  kinetics  includes  detailed  chemical  mechanism  involving  30 
species  {CsHg  ,0,  H,  O2,  N2,  H2,  CO,  OH,  H2O,  HO 2,  H2O2,  HCO,  CO 2,  CH,  CH2,  CH3,  CH4, 
C2H,  C2H2,  C2H3,  C2H4,  C2H5,  C2H6,  C3H5,  C3H6,  i-C3H7,  n-C3H7,  CH2O,  CH2OH,  CH3OH)  and 
75  elementary  chemical  reactions  taken  in  the  main  from  [5.4].  Only  rate  constants  of 
recombination  reactions  were  used  in  calculations.  If  in  the  reference  source  it  was  listed  rate 
constants  of  dissociation  reaction  so  rate  constant  of  recombination  reactions  is  calculated 
through  equilibrium  constant  of  this  reaction  and  it  is  approximated  by  Arrhenius  formula. 
Adiabatic,  fully  non-catalytic  model  of  gas-wall  thermal  and  chemical  interactions  was  used  in 
all  calculations. 

We  also  solved  the  system  of  chemical  equations  of  air  excited  by  non-selfmaintained 
discharge.  It  includes  equations  for  13  chemically  active  components  in  air,  ions  O^,  02  , 
negative  ions  O',  O2’,  O3',  atoms  O,  molecules  O2,  O3,  vibrationally  and  electronically  excited 
states,  including  electron  and  gas  temperatures  modified  in  chemical  reactions  as  in 

[5.2]. 

5.3  Results  of  calculations.  In  Fig.  5.2-3  one  can  see  calculations  of  propane  part  and  gas 
temperature  variation  dependencies  at  CsHg  burning  out  at  atomic  oxygen  addition  to 
stoichiometric  propane- air  mix.  The  calculations  were  made  for  initial  temperatures  of  T=  15  00 
and  2000  K  at  initial  atmospheric  pressure,  atomic  oxygen  part  in  the  mix  5=  0.0,  0.001,  0.01. 
Preliminary  calculations  showed  that  reasonable  temperature  range  of  propane  burning  out 
(ignition)  in  undisturbed  conditions  with  typical  times  of  t<200  ps  corresponds  to  initial 
temperatures  J>1800  K. 

The  presence  of  atomic  oxygen  in  the  mix  significantly  influences  propane  burning  out. 
Calculation  shows  that  appearance  of  0.1%  of  atomic  oxygen  in  the  mix  leads  to  ~1.5-2  times 
greater  reduction  of  propane  (ignition)  during  the  same  time  as  in  the  case  of  the  undisturbed 
mix.  Appearance  of  1%  of  atomic  oxygen  in  the  mix  leads  to  ~.5-6  times  greater  reduction  of 
propane  (ignition)  during  the  same  time  as  in  the  case  of  the  undisturbed  mix. 

In  Fig.5.4  a  -b  one  can  see  components  appearing  at  addition  of  atomic  oxygen  in 
reacting  gas.  Main  of  them  are  C3H6,  n-C3H7,  H2O.  At  initial  temperature  of  1500  K  and  part  of 
atomic  oxygen  5=  0.001  the  whole  process  of  molecular  transformation  takes  place  during  rch 
~40  ps,  and  it  decreases  to  Zch  ~20  ps  with  increasing  of  5  to  5=  0.01.  This  time  allows  to 
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estimate  the  plasma  area  width  needed  for  plasma  formation,  it  isZ^  ~  where  v/-is  flow 

velocity,  so  Lp~2  cm  at  vp-lGOO  m/s. 

These  calculations  show  that  if  the  temperature  in  the  mix  is  realized  by  external  plasma 
source  in  the  range  of  1800-2500  T,  and  the  dissociation  of  the  oxygen  molecules  leads  to  the 
appearance  of  0.01-01%  (of  the  total  mix)  oxygen  atoms  then  the  ignition  process  will  be 
realized  during  typical  times  smaller  than  the  gasdynamic  ones  (we  consider  typical  length  of 
plasma  region  created  by  non-self  maintained  discharge~0.02  m  and  M~2,  see  above). 

In  Fig.  5.5-7  one  can  see  temporary  dependencies  of  charged,  excited  and  active  particles 
created  in  the  plasma  of  non-selfmaintained  gas  discharge  excited  by  the  electron  beam  at  the 
excitation  rate  of  Q~10’’  cm'^  s'^  corresponding  a)  to  E-beam  power  of  about  ~0,6  W  and  a)),  at 
the  discharge  electric  field  of  10  ^  V/cm  and  b)  the  discharge  electric  field  of  10  V/cm. 

In  Fig.5.5  one  can  see  charged  particles  temporary  dependencies  in  such  a  plasma. 
Obtained  results  allow  to  estimate  current  density  in  realized  plasma  j-Qn^w^r,  where  e-electron 
charge,  Ue-electron  concentration,  w*-  electron  drift  velocity  (for  air  see  [5.7]),  obtained  results 
give 7~3  A/cm^  at  E=10  ^  V/cm,  and  600  A/cm^  at  E=10  V/cm. 

In  Fig.5.6-7  one  can  see  neutral  and  excited  states  temporary  dependencies  in  the 
plasmas. 

In  Fig.5.8  one  can  see  gas  and  electron  temperature  temporary  dependencies  in  the 
plasmas  created  by  non-selfmaintained  discharges. 

Calculations  show  that  favorable  for  combustion  enhancement  conditions  are  created  at 
Q~10*’  cm'^  s'*  and  the  discharge  electric  field  value  of  10  ^  V/cm  (electric  discharge  source 
power~10.0  kW).  They  correspond  to  E/N  value  of  about  10'*^  V  cm^ ,  it  is  near  the  threshold  of 
molecular  oxygen  dissociation.  High  temperature  realized  at  that  can  cause  the  ignition. 

The  calculations  also  show  high  number  (-10-20%)  of  vibrationally  excited  molecular 
oxygen  species  and  about  the  same  number  of  single  molecular  oxygen  (02d  =  02(A))  in  both 
cases.  These  components  can  reduce  additionally  the  ignition  time  and  reduce  necessary  power 
of  applied  plasma  source. 

Having  this  data  it  is  possible  to  estimate  the  power  W  of  the  discharge  section  of  2 
x5x5cm,  one  can  obtain  W~15  W  at  E=10  ^  V/cm,  and  W~30  kW.  The  last  value  can  be  easily 
decreased  to  W<10  kW  in  the  pulsed  mode  of  electric  field,  due  to  long  decay  of  plasmas  at 
temperature  T~2500K  [5.8]. 

Conclusions  to  Chapter  5. 

In  the  present  Report  we  considered  influence  of  excited  species  produced  by  the  non- 
selfmaintained  discharge  on  propane  ignition. 
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Calculation  shows  that  appearance  of  0.1%  of  atomic  oxygen  in  the  mix  leads  to  ~1.5-2 
times  greater  reduction  of  propane  (ignition)  during  the  same  time  as  in  the  case  of  the 
midistiirbed  mix.  Appearance  of  1%  of  atomic  oxygen  in  the  mix  leads  to  ~5-6  times  greater 
reduction  of  propane  (ignition  becomes  shorter  by  ~5-6)  during  the  same  time  as  in  the  case  of 
the  undisturbed  mix. 

Moderate  power  self-maintained  discharge  (<10.0  kW)  can  be  used  for  ignition  and 
ignition  enhancement  of  propane-air  stoichiometric  mix  due  to  generation  of  atomic  oxygen  and 
heating  of  the  mix  by  the  discharge.  High  amormt  of  vibrationally  and  electronically  excited 
states  can  considerably  reduce  necessary  power  value.  Application  of  discharges  combination: 
the  non  -selfmaintained  discharge  and  plasma  jet  (see  Fig.5.1)  at  low  non-selfinaintained 
discharge  power  (~20  W)  can  be  perspective  for  the  ignition,  since  vibrationally  and 
electronically  excited  oxygen  states  can  significantly  enlarge  reaction  rate  constants  and  can  be 
effective  at  interaction  of  incident  excited  mix  with  the  plasma  jet. 

More  precise  values  of  necessary  ignition  enhancement  power  value  require  detailed 
solution  of  non-selfinaintained  discharge  in  air-fuel  mix  including  rates  of  vibrationally  and 
electronically  excited  states  of  molecular  oxygen  and  hydrocarbon  molecules  reactions  and  the 
solution  of  gasdynamics  equations  in  the  definite  duct.  The  solution  of  the  first  problem  is 
complicated  because  of  necessary  rate  constants  absence.  The  solution  of  the  second  problem 
can  be  realized  in  approach  developed  in  this  Report  and  can  be  the  subject  of  further 
investigations. 


Fig.S.l  Schematics  of  the  channel  with  combination  of  discharges 

1  —  supersonic  M=2  nozzle,  2  —  collectors  of  preliminary  fuel  injection  (first  injector  set),  3  —  insulator  channel,  4  - 
collec  tors  of  main  fliel  injection  (second  injector  set),  5  —  1st  section  of  combustor,  6  -  2  ud  section  of  combustor,  7  - 
quartz  window,  8  -  third  fiiel  injector  set,  9— plasma  jet  generator,  10 -non-selfinaintained  discharge  generator 
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Fig.5.2a  Fig.5.2b 

Fig  5.2.  C3H8  burning  out  temperature  dependence  at  atomic  oxygen  addition  to  stoichiometric  propane-air 
mix. 

a)  initial  gas  temperature  T=1500  K,  b)  initial  gas  temperature  T=2000  K.  Part  of  oxygen  1-0.0,  2-0.001.,  3-0.01 


Fig.5.3a  Fig.5.3b 

Fig.5.  3.  Temperature  variation  at  C3H8  burning  out  temperature  dependence  at  atomic  oxygen  addition  to 
stoichiometric  propane-air  mix. 

a)  initial  gas  temperature  T=1500  K,  b)  initial  gas  temperature  T=2000  K.  Part  of  oxygen  1-0.0,  2-0.001.,  3- 

0.01 
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Fig.5.4a  Fig.5.4b 

Fig.  5.4.  Temporary  generation  of  active  particles  appearing  due  to  atomic  oxygen,  at  atomic  oxygen 
addition  to  stoichiometric  propane-air  mix.  Initial  gas  temperature  T=1500  K 

Part  of  atomic  oxygen  5  :  Fig4a  5=0.0,  Fig.4b  5=0.01 


Fig.5.5  Main  charged  particles  temporary  dependencies  in  the  plasmas  created  by  non-selfmaintained 
discharges.  Initial  gas  temperature  lOOOK,  Pressure  1  atm.  Upper  figure  -a)  -E-beam  source  power  -0.6  W, 
discharge  electric  field  of  -10^  V/cm;  lower  figure  -b)-E-beam  source  power  -0.6  W,  discharge  electric  field  of 
-10"^  V/cm  (p-positive  ions,  m-negative  ions) . 
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Fig.5.  6  Main  neutral  and  excited  states  temporary  dependencies  in  the  plasmas  created  by  non- 
selfmaintained  discharges.  Initial  gas  temperature  lOOOK,  Pressure  1  atm.  Upper  figure  -a)  -E-beam  source  power 
-0.6  W,  discharge  electric  field  of  -10^  V/cm;  lower  figure  -b)-E-beam  source  power  -0.6  W,  discharge  electric 
field  of -lO'*  V/cm.  (N2z-  N2(A),  N2b-  NjCB). 


76 


Fig.5.7  Main  vibrational  and  excited  states  temporary  dependencies  in  the  plasmas  created  by  non- 
selfmaintained  discharges.  Initial  gas  temperature  lOOOK,  Pressure  1  atm.  Upper  figure  -a)  E-beam  source  power 
-0.6  W,  discharge  source  power-1.0  kW;  lower  figure  -b)  E-beam  source  power  -0.6  W,  discharge  source  -10  kW. 
(02v,N2v-  effective  1-st  vibrational  levels  of  O2  and  N2,  molecules,  02d  -  02(A)). 
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Fig.5.8  Gas  and  Electron  temperature  temporary  dependencies  in  the  plasmas  created  by  non-selfmaintained 
discharges.  5.8a  -E-beam  source  power  -0.6  W,  discharge  source  power-1.0  kW;  5.8b  -E-beam  source  power  -0.6 
W,  discharge  source  -10  kW. 


Table  1.  Reactions  of  dissociation  -recombination  of  atomic  particles  [5.4] 

Rate  constant  form  =  AT ^  Qxp{—E  /  T)  ,  cmV(moP  s  ) 
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Table  2.  Rate  constants  of  exchange  reactions 

Rate  constant  form  =  AT ^  Qxp(—E  /  T)  ,  cmV(mol^  s  ) 
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H 

3.00el3 

0.0 

0. 

C2H4 

+ 

H02 

= 

C2H5 

+ 

02 

2.00el2 

0.0 

2513. 

C2H5 

+ 

H2 

= 

C2H6 

+ 

H 

5.40e02 

3.5 

2621. 

C2H5 

+ 

OH 

= 

C2H6 

+ 

0 

3.00e07 

2.0 

2573. 

C2H5 

+ 

H20 

= 

C2H6 

+ 

OH 

6.30e06 

2.0 

325. 

C3H5 

+ 

H2 

= 

C3H6 

+ 

H 

5.00el2 

0.0 

765. 

C3H6 

+ 

H02 

= 

C3H7I 

+ 

02 

1.00el2 

0.0 

2513. 

C3H7n 

+ 

H2 

= 

C3H8 

+ 

H 

1.30el4 

0.0 

4881. 

C3H7i 

+ 

H2 

= 

C3H8 

+ 

H 

1.00el4 

0.0 

4196. 

C3H7n 

+ 

OH 

= 

C3H8 

+ 

0 

3.00el3 

0.0 

2942. 

C3H7i 

+ 

OH 

= 

C3H8 

+ 

0 

2.60el3 

0.0 

2248. 

C3H7n 

+ 

H20 

= 

C3H8 

+ 

OH 

3.70el2 

0.0 

829. 

C3H7i 

+ 

H20 

= 

C3H8 

+ 

OH 

2.80el2 

0.0 

433. 

CH20 

+ 

H2 

= 

CH20H 

+ 

H 

3.00el3 

0.0 

0. 

CH20 

+ 

H02 

= 

CH20H 

+ 

02 

1.00el3 

0.0 

3619. 

CH20H 

+ 

H2 

= 

CH30H 

+ 

H 

4.00el3 

0.0 

3066. 

CH20H 

+ 

H20 

= 

CH30H 

+ 

OH 

1.00el3 

0.0 

853. 

C3H8 

+ 

02 

= 

C3H7n 

+ 

H02 

5.50e05 

1.5 

0. 
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C3H8  +  02 


C3H7i 


+  H02 


3.00e09 


0.5 


0. 


Table  3.  Reaction  rate  constants  in  non  selfmaintained  air  discharge,  T  in-°K,  [5.2], 


Number 

Reaction 

Rate  (cmVs,  cm^/s) 

/Cost  (eV) 

Comments 

En. 

relea 

se 

Energy 

defect, 

eV 

la 

C  +  O 2  — ^  O 2  C  C 

Ve=  141 

Konovalov 

lb 

c  +  O2  (v)  — ^  O2  c  C 

Ve=  141 

Konovalov 

Ic 

e  +  +  e  +  e 

Ve=  141 

Konovalov 

Idl 

C  +  N  2  — ^  N 2  ^  ^ 

Ve=43.1 

Konovalov 

ld2 

6  +  N 2  (^ )  — ^  ^2  ^  ^ 

Ve=43.1 

Konovalov 

ld3 

e  +  N^  )^N;+e  +  e 

Ve=  43.1 

Konovalov 

ld4 

e  +  N^(B^n^)^  N;  +e  +  e 

Ve=43.1 

Konovalov 

2al 

e  +  N2  N  +  N  +  e 

Vn  =33.6 

Konovalov 

N2 ,  Rydberg 

states,  with 

threshold 

12.25 

2a2 

e  +  )^N  +  N  +  e 

Vn*  =33.6 

Konovalov 

N2 ,  Rydberg 

states,  with 

threshold 

12.25 

2a3 

e  +  N^(B^n^)^N  +  N  +  e 

Vn*  =33.6 

Konovalov 

N2 ,  Rydberg , 

with  12.25  eV 

threshold 

2a4 

e  +  N2{v)  ^  N  +  N  +  e 

Vn*  =33.6 

Konovalov 

3a 

€  +  0 2  — ^  0  +  0  +  ^ 

Vo  =44.4 

Konovalov 

O2  all  from 

3b 

c  +  O2  (v)  — ^  0  +  0  +  ^ 

Vo  =44.4 

Konovalov 

O2  all  from 
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3c 

Vo  =44.4 

Konovalov 

O2  all  from 

3d 

6-\-02(b) — ^0+0+6 

Vo  =44.4 

Konovalov 

O2  all  from 

oyA^i:: 

4a  1 

^+Q  -^{^ffect)>t^ 

V  o2(a)=93.3 

Konovalov 

Oya'A), 

4a2 

e+Q(v) 

V  o2(a)=93.3 

Konovalov 

Oya^A), 

4a3 

^+Q  ~^(b)+e 

V  o2(b=386 

Konovalov 

0,ib% 

4a4 

e+C\(v)  -^(b)+e 

V  o2(a)=386 

Konovalov 

0,ib% 

4a5 

— ^C^fy^-hC 

V  O2(b=170 

Konovalov 

o,(^z:) 

4a6 

e+Q(v)—>Q{A)+e 

V  o2(A)=170 

Konovalov 

Sal 

e  +  02  -^0  +  0 

2-10'’-(300/Tef^ 

Smirnov 

gas 

+7.1 

5a2 

0  +  0  — ^  O2  +  € 

1.18-10'^'-  (T)^-^ 
exp(-80600/T) 

Gupta 

gas 

-7.1 

5bl 

O  +  M^O^+M  +  e 

2.0-10-+(T)‘-^-exp(- 

158000/T) 

Zamyshlyaev 

M=N2,  O2,  N, 

0,  NO 

gas 

-16 

5b2 

O^+M  +  e^O  +  M 

6  10'^^-(300/Te)‘  ^ 

Kossyi 

M=N2,  O2,  N, 

0,  NO 

gas 

+16 

5cl 

O  +  e^O^+e  +  e 

2.0-10-+(T)‘-^-exp(- 

158000/T) 

Zamyshlyaev 

Maikapar 

gas 

-16 

5c2 

O^+e  +  e^O  +  e 

9.05-10'^'’-(300/Tef^ 

Smirnov 

gas 

+16 
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6al 

e  +  NO^  ^N  +  O 

3.7-10'^-(300/Tef^^ 

Park 

gas 

+2.75 

6a2 

N  +  O^e^NO^ 

1.5-10'‘^-(T)-exp(-32  000/T) 

Park 

gas 

-2.75 

6b  1 

2.8-10-^-(300/Tef^ 

smimov 

gas 

+5.85 

6b2 

N  +  N^e  +  N^^ 

2.4-10-‘^-(T)'’'^-exp(-67 

300/T) 

Zamyshlyaev 

gas 

-5.85 

6c  1 

N  +  M^N^  +M  +  e 

2.0-10-^°-(T)‘'^-exp(- 

158000/T) 

M=N2,  O2,  N, 

0,  NO 

Zamyshlyaev 

gas 

-14.53 

6c2 

Kossyi 

gas 

+14.53 

6dl 

N^-e^N^  ^-e+e 

2.0-10-^°-(T)‘'^-exp(- 

158000/T) 

Zamyshlyaev 

gas 

-14.53 

6d2 

+e  +  e^N+e 

9.05-10'^°-(300/Tef^ 

Smirnov 

gas 

+14.53 

6fl 

N^+M  ^Nl+M  +  e 

2.0-10-^'’-(T)‘-^-exp(- 

158000/T) 

M=N2,  O2,  N, 

0,  NO 

Zamyshlyaev 

gas 

-15.6 

6f2 

N^^+M  +  e^N^+M 

6  10'^^-(300/Te)‘  ^ 

M=N2,  O2,  N, 

0,  NO  Kossyi 

gas 

+15.6 

6gl 

N 2  +  ^  ^  N 2  6  +  6 

2.0-10-^'’-(T)‘-^-exp(- 

158000/T) 

Zamyshlyaev 

gas 

-15.6 

6g2 

N 2  6 -\- 6  — ^  N 2  +  ^ 

9.05-10-^'’-(300/Tef^ 

Smirnov 

gas 

+15.6 

7a 

6  +  O2  +  O2  — ^  O2  +  O2 

^«=l-7-10-^‘- 

/  s 

Alexandrov 

gas 

-0.43 

7b 

O2  +  O2  — ^  O2  ^  0  2 

8.  10'*®  •exp(-6035/T) 

Mnatsakanya 

n,  Naidis 

gas 

+0.43 

8 

02(A)  +  O2  ^  e  +  O2  +  O2 

2.  lO'*" 

kossyi 

gas 

+0.55 

9 

02^  +  NO  ^  NO^  +  O2 

4.4  10'"^ 

kossyi 

gas 

+2.95 

10 

0  +  0 2  +  0 2  — ^  O3  +  0 2 

8.6  io'^'-r-‘^^ 

kossyi 

gas 

+1.04 

11 

0  +  O2  +  N 2  — ^  ^3  ^2 

5.6  lO'^^-T"" 

kossyi 

gas 

+1.04 

12a 

O2  {ci  A)  +  O3  — >  0  +  O2  +  O2 

9.7- 10’*^  exp(-1564/T) 

kossyi 

gas 

-0.06 

12b 

02(a'A;  +  A^^  A^O  +  O 

2.0  10"‘^  exp(-600/T) 

Kossyi 

gas 

+1.88 
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12c 

O2  ( ct  A  )  N  2  — ^  N  2  O2 

3.0  10'^' 

Kossyi 

gas 

+0.98 

12d 

O2  ( ct  A  )  O2  — ^  O2  ^2 

2.2  (T/300)"  '^ 

Kossyi 

gas 

+0.98 

12f 

O2  ( Cl  A)^0  — ^  O2  0 

7.0  10"‘^ 

Kossyi 

gas 

+0.98 

12g 

0^(a'A)  +  N0  ^O^+NO 

2.5  10'“ 

Kossyi 

gas 

+0.98 

13al 

N  2  0  2  — ^  TV  +  TV  +  0  2 

5.0- lO'^-  exp(-113200/T)- 
[1-  exp(-3354/T)] 

Krinonosova 

gas 

-9.76 

13a2 

TV  +  TV  +  O2  — ^  ^2  ^2 

8.27  exp(500/T) , 

T=1 00-600  K, 

1.9  10■^^  T=600-6300  K, 

Krinonosova 

gas 

+9.76 

13bl 

7V2  +  ^  ^  ^  +  ^2 

5.0- lO'^-  exp(-113200/T)- 
[1-  exp(-3354/T)] 

Krinonosova 

gas 

-9.76 

13b2 

+  A^2  ^  ^2  +  ^2 

8.27  lO'^^-  exp(500/T) , 

T=1 00-600  K, 

1.9  10■^^  T=600-6300  K, 

Krinonosova 

gas 

+9.76 

13cl 

5.0- lO'^-  exp(-113200/T)- 
[1-  exp(-3354/T)] 

Krinonosova 

gas 

-9.76 

13c2 

5.0  lO'^^-  exp(500/T) , 

T=1 00-600  K, 

1.14  10’^^  T=600-6300  K, 

Krinonosova 

gas 

+9.76 

13dl 

a^2  +  c>^a^  +  a^  +  c> 

l.MO'^-  exp(-113200/T)- 
[1-  exp(-3354/T)] 

Krinonosova 

gas 

-9.76 

13d2 

A^  +  A^  +  0^A^2  +  C> 

5.0  lO'^^-  exp(500/T) , 

T=1 00-600  K, 

1.14  10’^^  T=600-6300  K, 

Krinonosova 

gas 

+9.76 

13el 

1.8-10-^-  exp(-113200/T)- 
[1-  exp(-3354/T)], 

300-6000  K, 

1.4-10'^-  r^'^-exp(- 

Krinonosova 

gas 

-9.76 
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113200/T),  6000-20000  K, 

13e2 

N  +  N  +  N  ^  N^+  N 

1.8-10-^-  exp(- 113200/T)- 
[1-  exp(-3354/T)], 

300-6000  K, 

1.4-10-^-  r^-^-exp(- 
113200/T),  6000-20000  K 

Krinonosova 

gas 

+9.76 

13fl 

N^  +  0  ^  N  +  NO 

1.3  10'*®-  exp(-38000/T) 

Krinonosova 

gas 

-3.3 

13f2 

N  +  NO  ^  N^  +  0 

1.05-10-‘^-(T)^ 

200-4000  K 

Krinonosova 

gas 

+3.3 

14 

N  +  O  +  M^  NO  +  M 

1.8  io'^'-r“”^ 

Krivonosova 

gas 

+9.25 

15 

O2  +  0  +  O2  — ^  ^  +  O2 

3.0  lO'^*^ 

Smirnov 

gas 

+0.61 

15a 

0  2  +0  — ^  ^  +  0  2 

1.5  10'‘" 

Kossyi 

gas 

+0.61 

16 

0  2  +0  — >  0  2  0 

3.3  10'"^ 

kossyi 

gas 

+1.03 

17 

0  2  ^  3  — ^  ^3  ^  2 

4.  10"‘" 

Kossyi 

gas 

+1.57 

18 

0  +  O2  (A)  — ^  6  +  O3 

0^ 

0 

0 

Smirnov 

gas 

+0.5 

18a 

0  +  0  2  +  O2  — ^  0^  ^2 

1.1  10'^V300/T) 

Kossyi 

gas 

+0.54 

18b 

0  +  O2  — ^  ^2  ^ 

8.8  lO'*^  T'^"exp(- 
111837/T) 

Naidis 

gas 

-1.02 

18c 

0  +  O2  +  A^2  — ^  ^3  ^2 

1.1  10'^"(300/T) 

Kossyi 

gas 

+0.54 

18d 

0  +  O2  ^  y)  — ^  ^2  ^ 

1.  10'"^ 

Kossyi 

gas 

+2.01 

19 

O3  +  O2  — ^  C  +  O3  +  O2 

2.3  10'“ 

Samoilovich 

gas 

-2 

19a 

O3  +  0  — ^  O2  O2 

3.2  lO"'" 

Kossyi 

gas 

+0.56 

20 

0 2  +  O2  — ^  ^2  ^2 

4.2  lO’^  (300/T)‘'-^ 

Smirnov 

gas 

+11.77 

21 

0  2  0  2  ^  ^  ^2  ^2  ^ 

0.5- 10’^^  (300/T)'-^ 

gas 

+11.77 

21a 

O2  +  O2  ^2  — ^  ^2  ^2  0  -\-  0 

0.5- 10'^^  (300/T)'-^ 

gas 

+6.67 
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22 

1-  10'^(300/Tf^ 

gas 

+6.67 

23 

0~  +  0^  +M  ^  O  +  O^+M 

1.0- 10'^^  (300/T)'-^ 

gas 

+10.74 

24 

+  O 2  — ^  ^3  ^2 

2  10'^  (300/T)"  ^ 

gas 

+10.2 

25 

0^  +  O2  ”1”  ^  ^  ^3  ”1”  O2  4“  ^ 

1.0- 10’^^  (300/T)'-^ 

gas 

+10.2 

26 

O3  +  0^  — ^  0  +  0  +  O3 

2  10'^  (300/Tf  ^ 

gas 

+5.1 

27 

O^+O^+M^Oj^+O  +  O  +  M 

1.0- 10'^^  (300/T)'-^ 

gas 

+5.1 

27b 

0  +  0^  — ^  ^2  ^2 

2.0-1  O'* '  exp(-2275/T) 

Baulch 

gas 

+3.08 

28 

O"  +  NO +  0 

5-  10'\300/T)®-^ 

gas 

+7.79 

29 

0~  +N0^  +M^0  +  N0  +  M 

1.0-10'^^  (300/T)*-^ 

gas 

+7.79 

30 

Q-  +  TVO^  ^  iVO  +  C>2 

3-  10'\300/T)®-^ 

gas 

+8.82 

31 

01  +  NO^  +  M  ^  0^  +  NO  +  M 

1.0-10'^^  (300/T)*-^ 

gas 

+8.82 

32 

Q-  +  ^  O3  +  iVO 

2-  10'\300/T)®-^ 

gas 

+7.25 

33 

0~  +  NO^  +  M  ^  0^  +  NO  +  M 

1.0-10'^^  (300/T)'-^ 

gas 

+7.25 

34a 

6  -\-  6  -\-  O2  — ^  ^2  ^ 

9.05-10'^'*-(300/Tef^ 

gas 

+12.1 

34b 

^  ^  +  A^2  — ^  ^2  ^ 

9.05-10'^°-(300/Tef^ 

gas 

+15.6 

35 

e  +  e  + A^6>^  ->•  A^6>  +e 

9.M0'^°-(300/Tef^ 

gas 

+9.25 

36 

e+02^0  +0 

Reaction  is  approximated 

in  [2] 

under  E-field 

action 

below 

E/N=10‘‘^ 

V-cm^ 

Tes8.617-10' 

K. 

elect 

-1.47 

36a 

0  +0  — ^  ^  +  0  2 

5.0-10'*° 

kossyi 

gas 

+1.47 

37 

e+02  (ciA)^0  +0 

Reaction  constant  is 
approximated 

under  E-field 

action 

below 

E/N=10-‘’ 

V-cm" 

elect 

-2.45 
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Te=8.617-10- 

K. 

38 

e+02  (v—1)  -^0  +0 

Reaction  constant  is 
approximated 

under  E-field 

action 

below 

Te=8.617-10- 

5^'-p.  ^ 

K.  ’ 

elect 

-1.28 

39a 

e^02^  O2  (aA)+e 

Reaction  constant  is 
approximated 

below 

E/N=10'‘’ 

V-cm^ 

elect 

-0.98 

39aa 

^+(92(aA)^  O2 

Reaction  constant  is 
approximated 

elect 

+0.98 

39b 

e+OiCv)^  O2  (aA)+e 

K39 

elect 

40a 

e+02  -^0+  0  +e 

Reaction  constant  is 
approximated 

under  E-field 

action 

below 

E/N=10‘‘^ 

V-cm" 

Tes8.617-10' 

5^'-p.  Y' 

K. 

elect 

-5.11 

40b 

^+(92(v)  ^0+  0  ~^e 

K40 

■ 

40c 

e+02{A)  ^(9+  0  +6 

K40 

elect 

-4.12 

41 

e-^N2^N2(v) 

Reaction  constant  is 
approximated 

below  E/N= 
10-'^  V-cm^ 
Te=8.617-10' 

5^'-p.  ^ 

K.  ’ 

elect 

-.292 

41a 

e+N2(v)^N2  +e 

Reaction  constant  is 
approximated 

elect 

+0.292 

42 

e+02^02(v)  +e 

Reaction  constant  is 
approximated 

below 

E/N=10‘‘^ 

V-cm^ 

Te=8.617-10' 

K.  ’ 

elect 

-0.196 

42a 

e+02(v)^02  +e 

Reaction  constant  is 
approximated 

elect 

+0.196 

43 

N2(v)  +  M^N2+M 

7-10''*^-exp(-141/(T,K)'^^) 

(M-molec) 

gas 

+0.293 

44 

N2(v)  +  0^N2+0 

5-10''^-exp(-128/(T,K)'^^) 

gas 

+0.293 

45 

02(v)  +  ©2^  O2+  O2 

1.4M0'‘^- 

exp(-90/(T,  T<500, 

gas 

+0.196 
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mVs 

1.3M0'‘^ 

exp(-126/(T,  K)*^^),  m^/s 

T>500 

46 

C  +  ISl 2  — ^  ^2  (^)  ^ 

Ve=  14.5 

Konovalov 

(normalized 

on  pure  N2  in 

air) 

47 

e  +  O2  ^  O2  (v)  +  ? 

Ve=3.19eV 

(normalized 

on  pure  02  in 

air) 

48a 

^  +  O2  — ^  O2  +  2^ 

K48 

5-10''®<E/N 

V-cm" 

elee 

-12.1 

48b 

^  +  O2  (v)  — ^  O2  +  2^ 

K48 

elee 

-11.9 

48c 

e  +  O2  {^effec)  ^  ^2  + 

00 

elee 

-11.1 

48d 

^  +  A^2  — ^  ^2 

00 

elee 

-15.6 

48e 

6  +  A^2  (^)  — ^  ^2 

K48 

elee 

-15.3 

48f 

^  +  iV2(^'^:)^AA2"+2^ 

00 

elee 

-9.43 

48g 

e  +  N^(B^n^)^N^^+2e 

00 

elee 

-8.25 

49a 

e  +  N^^N^iA^i:l)  +  e 

Vn*  =258.0  eV 

Konovalov 

(normalized 

on  pure  N2  in 

air) 

Threshold 

E=6.17eV 

49b 

e  +  N^  -^N^(B^n^)  +  e 

Vn*  =249.0  eV 

Konovalov 

(normalized 

on  pure  N2  in 

air) 

Threshold 

E=7.35  eV 

50a 

e  +  N^^N^{A^K)  +  e 

Reaction  constant  is 
approximated 

3-10''®<E/N 

<E/N=2-10‘‘^ 

V-cm^ 

elee 

-6.17 

50aa 

e  +  N^(A^I^  )^e  +  N^ 

Reaction  constant  is 
approximated 

elee 

6.17 

50b 

e  +  N^  N^(B^n^)  +  e 

Reaction  constant  is 
approximated 

3-10''®<E/N 

<E/N=2-10‘‘^ 

V-cm^ 

elee 

-7.35 
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50bb 

e  +  N^(B^n^)^N^+e 

Reaction  constant  is 
approximated 

elee 

7.35 

51 

n,(a^e:)+n,(a^i:)^ 
N,(X'll)  +  N,(B'n^) 

2.6-10'‘® 

Kossyi 

gas 

+4.99 

52 

N,(A^Z:)  +  N,(A^i:)^ 
N,(X'll)AN,(B^nj 

l.MO-' 

N,(X'El) 

considered  as 

N2 

gas 

+4.99 

53a 

n,(a^i:)+n,(x'e;,v)^ 
N,(X'E;.w)  +  N,(B’n^) 

3.0-10'‘° 

N^(X^Il,x 

eonsidered  as 

N2 

gas 

-0.9 

53b 

N,(li:)+0,^N,+0+0 

2.54-10'‘^ 

Bauleh 

Kossyi 

gas 

+1.07 

53cl 

N,(A^Z:)  +  N, 

5.0-10-‘' 

Kossyi 

gas 

+1.18 

53c2 

N,(B^n^)^ 

N,(A^Z:)  +  hv 

1.5-10^  s'^ 

Kossyi 

gas 

+1.18 

53c3 

N^(B^nj  +  NO 
N,(A^Z:)  +  NO 

2.4- 10'*® 

Kossyi 

gas 

+1.18 

53c4 

N,(B^nj  +  o,^ 

N^+O  +  O 

3.0-10'‘® 

Kossyi 

gas 

+2.25 

54a 

n;  +  o,^o;  +  n. 

K54=2.86-10'''^  cmVs 

Charge 

exehange 

A^2  ion 

averaged 
over  N2  and 

Nlii 

appears  in 

N:^+2N2^> 

and 

Nl  +  O^^^ 

gas 

+3.4 

54b 

n;+o,{v)^o;+n. 

K54=2.86-10'''^  cmVs 

gas 

+3.4 

54c 

Nl+0,{K^J^0l+N, 

K54=2.86-10''‘’ cm^/s 

gas 

+2.4 

55al 

0 2  +  0 2  — ^  0  +  0  +  0  2 

3.7- lO’^-  exp(-59380/T)  • 

Krinonosova 

gas 

-5.1 
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[1-  exp(-2240/T)]  , 

1=300-4000; 

16.0-  (T)-^  -^exp(-59380/T) 

1=4000-20000; 

55a2 

0  +  0  +  O2  — ^  O2 

2.45  lO’^'- 

1=300-4000; 

Krinonosova 

gas 

+5.1 

55bl 

N  ^  0  +  0  +  N 

9.3- lO'^-  exp(-59380/T)  • 

[1-  exp(-2240/T)]  , 

T=3000-7000; 

Krinonosova 

gas 

-5.1 

55b2 

0  +  0  +  N  ^  0^+  N 

6.3  lO’^^- 

1=300-4000; 

Krinonosova 

gas 

+5.1 

55b3 

0^+N  ^  O  +  NO 

1.0-  10-'^ 

•  T-exp(-3150/T) 

Krinonosova 

gas 

-1.4 

55cl 

0^+  NO  ^  0  +  0  +  NO 

9.3- lO’^-  exp(-59380/T)  • 

[1-  exp(-2240/T)] 

T=300-4000; 

4.0-  (T)'^  ^exp(-59380/T) 

T=4000-20000; 

Krinonosova 

gas 

-5.1 

55c2 

0  +  0  +  NO  ^  0^+  NO 

6.3  lO’^^- 

1=300-4000; 

Krinonosova 

gas 

+5.1 

55dl 

0  2  +  0  — ^  0  +  0  +  0 

1.3-10-^-  exp(-59380/T)  • 

[1-  exp(-2240/T)] 

T=300-4000; 

5.8  10-^(T)-^-^exp(- 
59380/T) 

T=4000-20000; 

Krinonosova 

gas 

-5.1 

55d2 

0  +  0  +  0  — ^  0  2  +  0 

8.8  lO'^'- 

T=300-4000; 

Krinonosova 

gas 

+5.1 

55el 

0  2  +  N  2  — ^  0  +  0  +  N  2 

9.3- lO’^-  exp(-59380/T)  • 

[1-  exp(-2240/T)] 

Krinonosova 

gas 

-5.1 

55e2 

0  +  0  +  7^2  — ^  ^2  ^2 

2.76  lO'^^-  exp(720/T) 

Krinonosova 

gas 

+5.1 
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56al 

NO  +  0  2  — ^  TV  +  0  +  O2 

8.7  lO-'" 

•  exp(-76  000/T) 

Krinonosova 

gas 

-6.5 

56a2 

N  +  0  +  O2  — ^  NO  +  O2 

2.8  (T)'°-^ 

Zamyshlyaev 

gas 

+6.5 

56bl 

NO +  N^^N  +  0  +  N^ 

8.7  10"^ 

•  exp(-76  000/T) 

Krinonosova 

gas 

-6.5 

56b2 

N  +  O  +  N^^  NO +  N^ 

1.8  lO’^'- 

Krinonosova 

Sovpadaet  s 

14 

gas 

+6.5 

56cl 

NO +  0^N  +  0  +  0 

1.7  lO'^ 

•  exp(-76  000/T) 

Krinonosova 

gas 

-6.5 

56c2 

N  +  0  +  0^  NO +  0 

5.6  lO'^’-  (T)'°-^ 

Zamyshlyaev 

gas 

+6.5 

56c3 

NO  +  0  ^  N  +  0^ 

2.5  10'‘^-  (T) 

•  exp(-19  500/T) 

Krinonosova 

gas 

+1.4 

56dl 

NO +  N^N  +  0  +  N 

1.3  10"^ 

•  (T)'  ‘••^exp(-75530/T) 

Zamyshlyaev 

gas 

-6.5 

56d2 

N  +  O  +  N^  NO +  N 

5.6  lO'^’-  (T)-®-^ 

Zamyshlyaev 

gas 

+6.5 

56d3 

NO  +N  ^0  +  N^ 

1.05  lO'*^ 

•  (Tf-^ 

T=200-4000  K 

Krinonosova 

gas 

+3.26 

56el 

NO  +  NO  ^  N  +  0  +  NO 

1.7  10"^ 

•  exp(-76  000/T) 

Krinonosova 

gas 

-6.5 

56e2 

N  +  0  +  NO  ^  NO  +  NO 

5.6  lO'^’-  (T)-®-^ 

Zamyshlyaev 

gas 

+6.5 

58 

N  +  NO  ^  N^+0 

1.0  lO’*^-  (T)®  -^ 

Krinonosova 

gas 

+3.28 

62 

N  +  0^^  NO  +0 

4.1  10’*^  •exp(-3200/T) 

Bauleh 

gas 

+1.4 

63 

N  +  0,  ^  NO  +  0^ 

5.5  10'*^ 

Bauleh 

gas 

-3.3 

64 

N,(J^rj^O^N,  +0 

2.1  10'“ 

Kossyi 

gas 

+6.17 

65 

N^(jirj+O^NC^N 

7.0  10'*^ 

Kossyi 

gas 

+2.91 

66 

N2 

2.5  10-‘'’-exp(-50391/T) 

krivonosova 

gas 

-3.9 

67 

0-  +NO^e  +  NO^ 

2.6  10'*" 

Kossyi 

gas 

+1.64 

68 

O2  (b )  ^  0^  — ^  0  +  O2  +  O2 

1.8-10'‘' 

kossyi 

gas 

+0.6 

69 

0^(b)  +  N^  N^+0^(a^A) 

4.9  10-‘^-exp(-253/T) 

Kossyi 

gas 

+0.66 
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71 

O2  (b )  -\-  0 2  — ^  O2  "b  O2  ( ^ 

4.3  exp(-241/T) 

Kossyi 

gas 

+0.66 

72 

0yb)  +  0^0ya'A)  +  0 

8.0 

Kossyi 

gas 

+0.66 

73 

Oyb)  +  NO^  Oya'A )  +  NO 

4.0  10'*^ 

Kossyi 

gas 

+0.66 

74 

O2  (  A)  -\-  0 2  — ^  O2  (b )  F  0 2  (b ) 

2.9  10"'^ 

Kossyi 

gas 

+1.22 

75 

oyA)  +  o^oyb)  +  o(D) 

9.0  10"'^ 

Kossyi 

gas 

+0.89 

76 

OyA)  +  N^  -^N^+Oyb) 

3.0  lO'*^- 

Kossyi 

gas 

+2.86 
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Chapter  6.  Stationary  and  pulsed  plasma  jets  for  ignition 
6.1.  Technical  starting  of  the  plasma  generator  with  the  divergent  conical  nozzle 
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In  accordance  with  the  technical  task  and  the  calendar  plan  of  works  technical  starting  of 
the  plasma  generator  with  the  divergent  conical  nozzle  has  been  realized.  (Main  elements  of  the 
plasma  generator  constinction  were  considered  in  the  Report  [1]). 

Nitrogen  and  argon  were  used  plasma  fomiing  gases.  Modes  of  the  plasma  generator 
work  have  been  studied,  its  Ampere-Volt  characteristics,  distributions  of  plasma  jet  parameters 
have  been  studied  with  respect  to  the  sort  of  the  gas  and  its  flow  rate.  The  gas  flow  rate  varied  in 
the  range  1  —  1,5  g/s.  Strong  erosion  of  the  electrode  surface  arises  at  gas  flow  rates  smaller  than 
1  g/s.  Main  series  of  the  experiments  was  carried  out  at  gas  flow  rate  of  1,25  g/s.  This 
corresponds  to  nitrogen  flow  rate  of  ~  75  liter/minute  and  argon  flow  rate  of  -  40  liter/minute. 


Starting  block 


Power  source 


Exhaust  (forced) 
ventilation 


Workiup 


Plasma 


Basement 


Fig.l 


hi  Fig.6. 1  one  can  see  a  photo  of  the  set  up.  Main  elements  of  the  set  up  are  the  follows: 
1 .  Power  source  with  power  <15  kW; 
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2.  Starting  block,  it  fonns  an  ignition  pulse  at  amplitude  of  ~  5  kV,  duration  of  an  ignition 
pulse  is  about  of  100  ps; 

3.  Basement,  in  which  placed  systems  of  water  cooling,  current  stabilizing,  gas  delivery  and 
consumption; 

4.  Working  chamber  with  the  plasma  generator  ( the  photo  is  made  without  of  protecting 
glass  —  a  filter) 

5 .  System  of  forced  ventilation. 


Working 

chamber 

Nozzle 

Plasma 


Water  cooling 
system  element 


The  plasma  generator  appearance  in  Fig-6-2  r  chamber  is  represented  in  Fig.6.2. 

hi  Fig.  6.3  -  Fig.  6.5  one  can  see  photos  of  the  plasma  jet  at  different  values  of  the 
inputted  power  (Fig.  6.3  -  3  kW;  Fig.  6.4  -  4,5  kW;  Fig.  6.5—  6,5  kW) 


Fig.  6.3 


Fig.6.4 
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Fig.6.5 

First  two  photos  were  made  with  the  application  of  dense  neutral  light  fdter.  Non 
symmetrical  picture  of  the  jet  luminescence  can  be  described  by  large  vortex  component  (see 
part  2  of  the  Report  about  the  vortices  in  the  channel)  of  the  velocity.  (Remember  that  in  the 
present  plasma  generator  design  the  tangential  swirling  is  used  for  the  flow  stabilizing  and  for 
the  gasdynamic  protection  of  electrodes.  At  that  the  vortical  velocity  component  Vr  equals  to  the 
axial  component  of  the  flow  velocity  Va,  and  the  angle  between  them  is  a  «  45°). 

6.2.  Determination  of  the  gas  temperature  in  the  discharge  created  by  the  plasma  generator 

with  divergent  conical  nozzle 

Information  about  the  gas  temperature  Tg  is  important  for  an  analysis  of  plasma 
processes.  The  most  simple  situation  in  case  of  molecular  plasma  corresponds  to  local 
thermodynamic  equilibrium  when  Tg  =  Ty  =  Tr,  here  Ty  and  Tr  -  are  the  vibrational  and 
rotational  temperatures,  respectively.  Ty  and  Tr  -  are  the  parameters  determining  vibrational 
and  rotational  level  populations  in  the  Boltzmann  function  characterizing  molecular 
distribution  function.  Both  Ty  and  Tr  can  be  obtained  from  radiation  spectra  of  the  plasma. 

Tg  in  non-equilibrium  plasmas  also  can  be  determined  on  a  basis  of  Tr  measurements  [2]. 
Plasma  radiation  spectra  give  information  on  rotation  temperatures  T*r  of  molecules  in  excited 
electronic  states.  T*r  dependence  of  the  spectral  components  hk,  belonging  to  the  rotational 
structure  of  spectral  lines  corresponds  to  the  transition  from  the  1-st  and  the  k-th  states  is 
described  by  the  formulae 

(ke  /B*)ln(Iik/i  v^k)  =  C  -  j(j+l)/T*R, 

here  ke  is  the  Boltzmann  constant,  j  -  full  molecular  angular  momentum,  i  -  quantum 
mechanical  intensity  coefficient,  nik  -  radiation  frequency  of  the  quantum,  C  -  constant  One 
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can  determine  T*r  from  the  graph  of  the  function  (ke  /B*)ln(Iik/i  v'^ik)  through  j(j+l)  and 
measurement  of  its  slope  tg  a: 

T*r  =  1/tg  a. 

The  probability  of  rotational  energy  changing  is  negligibly  small  for  a  stable  molecule, 
and  rotational  distribution  functions  over  the  excited  and  ground  states  of  molecules  are 
identical,  Tr  =  T*r.  Intense  exchange  of  translational  and  rotational  energies  at  collisions  leads 
to  equality  of  Tg  and  Tr  of  molecules  in  the  ground  electronic  state.  But  there  are  some 
features  in  real  conditions,  which  have  to  be  accounted.  If  molecules  take  part  in  plasma 
chemical  reactions  then  some  part  of  their  activation  energy  can  transfer  to  the  rotational 
energy,  this  leads  to  the  distortion  of  the  initial  rotational  distribution  function,  Tr  ^  T*r  in  this 
case.  This  is  typical  for  chemically  active  radicals.  Another  source  of  errors  is  connected  with 
the  fact  that  one  determines  not  Tr  from  spectra  but  B*/Tr,  here  B*  is  the  rotational  constant 
of  the  upper  state  at  the  radiation  transition.  The  rotational  energy  distribution  can  stay  the 
same  at  excitation  of  this  state  by  the  electron  collision  but  it  corresponds  to  the  temperature 

Tr  =  Tg  B*/B°, 

here  B°  is  the  rotational  constant  for  the  ground  state.  So  it  is  necessary  to  choose  stable 
molecules  with  B*  »  B°  for  determining  of  Tg. 

If  one  uses  the  second  positive  nitrogen  system  then  in  this  case  one  can  suppose  vV  to 
be  constant  (the  corresponding  error  is  ~  0,5  %,  and  it  can  be  neglected),  B*  =  1,826;  so  the 
formulae  for  Tr  has  the  following  form  in  this  case: 

0,89  ln(Iik/i)  =  Const  -  j(j+l)/TR. 


A  set  up  for  spectral  measurements  of  plasma  parameters  is  represented  in  Fig.8.1.6.  The 
discharge  radiation  was  focused  to  the  inlet  slot  of  the  spectrograph  STE  -1  with  a  help  of  quartz 
lens  (1)  with  the  focal  distance  f  =  50  cm  with  scaling  down  1  :  2.  The  device  dispersion  was~ 
0,35  nm/mm  in  the  range  of  waves  300  -  400  nm.  The  discharge  image  was  oriented  along  the 
inlet  spectrograph  hole  for  obtaining  of  spatially  resolved  radiation  spectrum  . 
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discharge 


lens 


Fia.  6.6.  A  scheme  of  radiation  snectra  detection. 

For  example  in  Fig.8.1.7  is  represented  the  discharge  radiation  spectrum  in  the  range  of 
wavelengths  7.  =  360  -  400  nm  for  a  discharge  with  electric  inputting  power  3  kW.  In  this  case 
the  electrode  erosion  was  rather  weak.  It  is  probable  that  CN  molecules  in  present  experiments 
are  generated  at  recombination  of  nitrogen  and  carbon  atoms.  Atomic  carbon  is  the  product  of 
interaction  of  the  plasma  with  the  inter  electrode  partitions  that  contain  organic  elements. 

At  low  inputting  power  it  is  possible  to  detect  the  gas  temperature  (nitrogen  in  this  case) 
on  a  basis  of  non  resolved  CN  bands  with  the  threshold  wavelengths  3883  A°  and  3871  A°  [3]. 
These  methods  are  based  on:  1)  ratio  of  radiation  integrals  of  two  bands;  2)  ratio  of  the 
threshold  intensities  of  these  bands;  3)  rotation  transition  intensities  distribution  of  the  band 
(0;0)  with  the  quantum  wavelength  3883  A°. 

The  spectrum  changes  at  increasing  of  the  inputted  power:  intense  lines  of  copper  and 
manganese  appear.  The  temperature  in  this  case  can  be  determined  on  a  base  of  relative 
intensities  of  chosen  line  pairs. 


rtPii; 
Ha  Cv  'v") 
A.,  nm: 
CN  (V  'V") 

nm: 

Nj'  Cv 


3H7.1  388-3 
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Fig.6.7.  Nitrogen  plasma  jet  radiation  spectrum  at  the  inletting  power  3  kW. 


Analysis  of  radiation  spectra  has  shown  that  the  gas  temperature  of  the  plasma  jets  in  the 
investigated  range  of  inputted  power  in  the  plasma  relaxation  region  (behind  the  nozzle  cut) 
lies  in  the  range  Tg  =  (3500  -  6600°)  K.  Its  longitudinal  dependence  is  not  monotonous.  Axial 
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temperature  distributions  show  that  the  gas  temperatiue  rises  with  the  cunent  rise  and  it 
decreases  down  the  flow  (see  Fig. 6. 8). 


Tg,K 


Fig.6.8.  Longitudinal  distribution  of  the  gas  temperature  of  the 
nitrogen  jet:  Inputted  power  1  -3  kW;  2  —  4,5  kW,  3-  6,5  kW. 

Spatial  distribution  of  electron  concentration  can  be  obtained  with  good  accuracy  by  the 
Saha  equation  (see  [4])  at  known  gas  temperatiue  distribution.  In  such  a  way  the  electron 
concentration  near  the  nozzle  cut  is  several  units  10  cm'  in  the  case  of  the  inputted  power  6,5 
kW  and  it  is  in  agreement  with  results  of  our  calculations. 

6.3.  Results  of  experiments  carried  out  with  the  plasma  generator  with  divergent  conical 

nozzle 

hi  accordance  with  the  technical  task  and  the  calendar  plan  of  works  investigations  of 
main  parameters  of  the  plasma  flowing  out  of  the  stationary  plasma  generator  with  the  divergent 
conical  nozzle  have  been  made. 

Modes  of  the  plasma  generator  operation  in  the  wide  range  of  the  initial  conditions  have 
been  studied. 

The  mass  flow  rate  of  a  gas  was  varied  in  the  range  0,7  —  2  g/s,  the  power  characteristics 
were  varied  in  the  range  3-7  kW. 

It  was  shown  that  the  gas  temperature  at  the  outlet  from  the  nozzle  changes  in  the  range 
from  3500  K  to  6000  K  at  changing  of  the  inputted  power  in  the  range  from  3  to  6,5  kW. 

The  typical  values  of  the  electron  density  lie  in  the  range  7-10^^  -  5*10^®  cm  ^  as  it  was 
obtained  by  the  scanning  electric  probe  fransversally  to  the  axis  in  the  mode  of  the  eontinuous 
medium. 

It  was  shown  that  the  strong  erosion  of  the  electiode  surface  takes  plaee  at  the  mass  flow 
rates  smaller  than  1  g/s  even  at  low  values  of  the  inputted  power  (P<  4  kW),  it  leads  to  their  rapid 
destruction. 
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The  detailed  analysis  of  the  planned  experiments  conditions  (7-10  quarters)  on  the 
stationary  supersonic  gaseous  and  plasma  jets  shows  that: 

1.  developed  and  put  into  the  exploitation  powerful  stationary  plasma  generator  is  highly 
effective  plasma  source,  which  stably  works  at  high  values  of  the  power  put  in  (up  to  10  kW)  and 
at  high  mass  flow  rates  of  the  plasma  forming  gas  (in  particular  nitrogen  and  argon).  It  can  be 
successfully  applied  for  the  plasma  chemical  activation  of  the  high  speed  (supersonic)  flows  of 
the  air-fuel  mixture  in  the  large  scaled  stationary  engine  devices,  which  transversal  sizes  are  no 
less  than  20  -  30  cm.  Namely  this  type  of  the  plasma  generator  can  be  recommended  for  the 
applications  in  the  set  ups  of  these  devices; 

2.  At  the  application  of  the  engine  model  with  the  typical  size  of  2-3  cm  (developed  and 
manufactured  in  MSU)  and  working  in  the  quasi-stationary  mode  (the  duration  of  the  working 
cycle  is  ~2-3  s  )  the  using  of  such  stationary  high  temperature  source  is  unlikely  advisable  since: 

a)  the  high  thermal  loads  can  lead  to  the  rapid  destruction  of  the  supersonic  channel  walls; 

b)  it  will  substantially  change  the  gasdynamic  structure  of  the  main  supersonic  flow  due  to 
large  additional  injection  of  the  plasma  forming  gas; 

c)  it  will  lead  to  the  uncontrolled  experimental  conditions  over  the  initial  conditions  (rapid 
uncontrolled  pressure  increase  in  the  receiver) 

d)  it  is  connected  with  large  difficulties  at  the  synchronizing  of  the  main  supersonic  tract 
work  ,  the  plasma  source  and  the  delivery  device  of  the  plasma  forming  gas  (rather  long  time, 
about  of  10  -  20  s,  is  required  for  the  plasma  generator  to  reach  the  stable  mode  of  the 
operation). 

6.4.  Small  sized  quasi  stationary  plasma  generator 

In  connection  with  conclusions  of  #8.1.1.  the  small  sized  quasi  stationary  plasma  generator 
has  been  developed  by  us  and  put  to  the  manufacturing.  The  nozzle  geometry  corresponds  to 
those  of  the  main  already  developed  construction.  The  typical  parameters  of  this  plasma 
generator  are: 

a.  Sizes  -  the  external  diameter  is  16  mm,  the  length  is  80  mm; 

b.  the  voltage  pulse  duration  -  5  -  6  ms; 

c.  -280  V; 

d.  the  discharge  current  maximum  value  120  A; 

e.  the  discharge  current  pulse  form  is  the  bell  type 

The  plasma  generator  general  scheme  is  represented  in  Fig.6.9  . 
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Fig.6.9  -  the  nozzle,  2  -  the  body,  3  —  the  dielectiic,  4  —  the  joint,  5  —  the  centia] 
electrode,  5-  the  coaxial  cable. 


The  plasma  generator  represents  the  cylindrical  axially  symmetric  construction.  Its  main 
elements  are  made  of  copper  (  the  nozzle,  the  body,  the  mormting  joint  for  the  plasma  generator 
connection  with  the  coaxial  cable  ).  The  central  (axial)  electrode  is  made  of  the  tungsten.  The 
organic  glass  or  the  heat  proof  ceramics  are  used  for  the  dielectric.  The  thread  Ml  6  is  foreseen  in 
the  front  part  of  the  plasma  generator  nozzle.  It  is  required  for  the  adjoining  of  the  plasma 
generator  to  the  supersonic  tract  or  to  the  intermediate  receiver  chamber. 

The  plasma  generator  is  placed  in  the  dhect  vicinity  to  the  supersonic  flow  on  the 
twisting  device,  which  allows  to  investigate  the  interaction  of  the  high  speed  gaseous  and  plasma 
flows  directed  at  different  angles. 

He  plasma  generator  is  connected  with  the  power  soruce  with  the  outlet  voltage  of  a  300 
V.  The  voltage  pulse  duration  can  be  varied  in  the  range  from  1  to  10  ms.  The  energy  stored  in 
the  capacity  storage  is  about  of  145  J,  and  it  corresponds  to  the  power  put  of  145  kW  to  the  load 
at  the  pulse  duration  of  «  1  ms. 

The  generator  starling  is  realized  by  the  voltage  high  voltage  pulse  with  the  amplitude  of 
30  kV.  Tire  possibility  of  the  low  voltage  synchronizer  prrlse  application,  which  comes  from  the 
starting  scheme,  allows  to  start  the  plasma  generator  in  any  time  period  with  respect  to  the 
general  pneumovalve  starting  the  supersonic  charmel. 

6.5.  Plasma  parameter  measurements  of  small  sized  quasistationary  plasma  generator 

For  generation  of  quasi-stationary  plasma  jets  the  small  sized  energy  storage  with  total 
capacity  of  C  =3200  pF  has  been  applied.  The  design  of  the  energy  storage  controlling  block 
allowed  to  work  at  different  fixed  values  of  the  stored  energy,  namely,  50,  100  and  200  J.  The 
operating  voltage  was  ~  300  V.  The  ciment  pulse  dmation  was  2,  12  and  14  ms.  The  maximum 
discharge  current  value  was  150  A,  the  ciuient  pulse  form  was  the  bell  type. 
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The  general  plasma  generator  design  was  diseussed  in  details  in  the  previous  part.  The 
plasma  generator  represented  the  cylindrieal  axisymmetrical  eonstruction,  its  main  parts  were 
made  of  eopper  (the  nozzle,  the  body,  fixing  joint  for  conneetion  of  the  plasma  generator  and  the 
eoaxial  eable).  The  eentral  axial  electrode  is  made  tungsten.  Organic  glass  was  used  as  the 
dielectric  material.  Starting  of  the  scheme  and  detection  devices  was  realized  by 

synchronizing-pulse  coming  from  the  generator  G5-15. 

Investigation  of  the  plasma  injecting  to  air  was  undertook  at  the  invariable  plasma 
generator  working  chamber  geometry.  Only  imputed  energy  to  the  discharge  was  varied.  Spectra 
detection  was  realized  by  the  spectrograph  STE  -1. 

Schlieren  method  was  applied  for  plasma  flow  character  (laminar  or  turbulent)  determination. 
Shadow  pictures  were  obtained  with  a  help  of  the  super  high-speed  photo  recorder. 

In  Fig.6.10  one  can  see  the  integral  photo  of  the  plasma  jet  in  the 
motionless  air  at  the  pressure  of  1  atm.  This  plasma  jet  has  almost 
right  cylindrical  form,  it  is  surrounded  by  weakly  luminescent  low 
temperature  “coat”  (but  the  area  directly  adjoining  the  jet,  where  the 
plasma  temperature  is  the  highest).  The  plasma  boundaries  are  not 
disturbed,  the  flow  is  close  to  the  laminar  one.  Investigations  of  the  jet 
dynamics  have  shown  that  the  plasma  flowing  out  from  the  injector 
nozzle  takes  place  with  considerable  delay  time  (about  of  800  -  900 
ps)  with  respect  to  the  current  pulse  beginning.  The  discharge  was 
developing  inside  the  plasma  generator  working  chamber  during  this 
time  period. 

Fig.  6.10.  Plasma  jet 

The  plasma  flow  character  was  investigated  with  a  help  of  the  shadow  method.  In  Fig.6.11  one 
can  see  fragments  of  the  frame-accurate  shadow  image  corresponding  to  the  initial  (a,  b  -  rising 
of  the  discharge  current)  and  the  final  (about  5  ms)  phases  of  the  process.  They  were  obtained  at 
the  stored  energy  value  of  50  J. 
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Fig.6.1 1.  Fragments  of  the  jet  shadow  image 
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One  can  see  that  the  flow  is  the  laminar  at  the  mitial  stage  of  the  process.  The  noticeable 
tnrbnlization  of  the  plasma  smioundiug  “coat”  is  obsei'ved  in  the  end  of  the  energy  put  pulse. 

Plasma  parameters  diagnostics  was  carried  out  by  the  optical  methods  over  the  jet 
radiation  spectrum.  It  proved  out  to  be  that  with  respect  to  the  power  input  there  were  the  lines  of 
CN  and  carbon  C2  in  different  ratios,  the  continuous  spectrum  detected  from  regions  near  the 
plasma  generator  nozzle,  lines  of  the  Ha,  Hp,  Hy,  copper,  sodium,  carbon,  etc.  The  hydrogen 
lines  Ha,  Hp,  Hy  were  strongly  spread,  at  that  the  spreading  value  depended  on  the  spatial 
coordinate  (in  the  longitudinal  and  the  transversal  directions)  of  the  jet. 

The  plasma  temperatme  was  measmed  either  over  the  lines  of  CN  with  the  quantum 
wavelength  X,  =  383,3  nm  (the  line  0;0)  or  over  copper  lines  with  the  wavelength  X  =  510,55  nm, 
5 15,32  nm  and  521 ,82  mn. 

The  sample  of  the  specfrum  image  of  the  line  (0;0)  is  represented  in  Fig.  6. 12. 


The  electron  concentrations  were  deteimined  by  Stark  spreading  of  the  hydrogen  line  Hq 
and  by  the  electric  probe  method. 

hi  Fig.6.13  and  6.14  one  can  see  the  dependencies  of  the  maximum  values  of  the  plasma 
temperature  and  the  electron  concentrations  near  the  plasma  generator  nozzle  cut  with  respect  to 
the  energy  stored  in  the  power  storage.  One  can  see  that  the  plasma  temperature  rises  from  ~  4 
kK  to  ~  6,5  kK  in  the  investigated  ranges  of  the  stored  energy,  at  the  same  time  the  electron 
concentration  rises  from  ~  2-10  cm"  to  5,5-10  cm"  at  the  energy  W  variation  from  50  J  to 


200  J. 
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Fig,  6,13 


Fig..6.14 


,6.15 

for  two  values  of  energy  stored  in  the  power  storage:  100  J  (the  curve  2)  and  200  J  (the  curve  1). 
One  can  understand  fiom  the  represented  data  that  the  temperature  monotonously  decreases  with 
increase  of  the  distance  from  the  plasma  generator  cut,  and  at  a  distance  of  about  10  cm  it  is  ~3 
kK  for  100  J  and  it  is  5,5  kK  at  200  J. 


Fig,6.15 


The  stark  spreading  of  the  hydrogen  lines  is  much  greater  than  the  spreading  due  to  other 
factors  at  the  electr  on  concentration  Ne>  10^^  cm'^.  Spreading  of  spectral  lines  in  conditions  of 
om'  experiments  can  take  place  due  to  the  Doppler  effect  (~  0,01  mn),  external  Stark  effect  (~ 
0,01  nm),  spectrograph  instrumental  function  (~  0,01  nm),  but  the  spreading  of  the  line  Ha 
detected  in  the  experiments  varied  from  0,1  to  3,5  nm.  The  choice  of  the  line  Ho-  for  the  electron 
concentration  determination  was  conditioned  by  the  following:  the  accuracy  of  the  line  spreading 
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theoretical  calculations  is  sufficiently  high  (~  10%)  and  the  given  line  corresponds  to  the 
transition  between  considerably  low-lying  levels.  So  it  is  easily  excited  m  plasmas  even  the 
hydrogen  is  present  as  small  admixture.  In  Fig. 6. 16  one  can  see  the  axial  distribution  of  the 
election  density  for  two  values  of  energy  stored  in  the  power  storage:  100  J  (the  cmve  2)  and 
200  J  (the  curve  1). 


One  can  see  that  the  electron  concentration  at  distances  in  the  range-l  cm  decreases 
comparably  slowly,  and  then,  starting  fiom  distances  about  1-  1,5  cm  it  decreases  practically 
linearly.  If  one  considers  that  at  larger  distances  -5  —  10  cm  it  decreases  linearly  with  the 
longitudinal  coordinate  then  he  obtains  the  electron  concentration  —  10^^  cm'^  at  the  distance  of 
10  cm.  This  was  confirmed  by  results  of  the  probe  measurements. 

From  the  analyses  of  the  experimental  data  one  can  make  the  following  conclusions: 

The  developed  and  applied  generator  of  the  quasi-stationaiy  plasma  jets  reprodirces  main 
plasma  characteristics  of  the  powerful  plasma  generator  with  the  expendable  power  up  to  15 
kW. 

The  spatial  extension  of  the  plasma  jet  of  the  quasi-stationary  plasma  generator  reaches  10 
cm. 

Conclusions  to  chapter  6 

Analysis  of  the  powerful  plasma  generator  application  possibilities  for  the  investigation 
of  the  plasma  and  gas  crossflows  have  been  made  on  a  basis  of  the  obtained  results. 

It  was  established  that  such  an  experiment  can  be  undertaken  at  the  large  sized  supersonic 
devices  with  the  cross  section  no  less  than  20  cm. 
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It  was  shown  that  the  developed  plasma  generator  is  the  highly  effective  plasma  source 
stably  working  at  high  energy  put  in  and  large  mass  flow  rates  of  the  plasma  forming  gas.  It  most 
probably  can  be  successfully  applied  at  the  plasma  chemical  activation  of  the  high  speed 
(supersonic)  flows  of  the  fuel-air  mixtures. 

For  the  realization  of  the  experiment  in  the  small  sized  devices  it  is  advisable  to  use  the  small 
sized  quasi  stationary  plasma  generators  with  the  external  controlled  starting.  Such  plasma 
source  has  been  developed  and  manufactured. 

The  range  of  plasma  parameters  of  the  quasi-stationary  plasma  generator  allows  to  assert  that 
this  type  of  the  plasma  source  is  of  use  for  ignition  of  the  fuel  supersonic  flows  both  in  the  region 
of  the  flow  stabilization  and  in  the  region  of  the  main  fuel  flow. 

This  plasma  source  has  applicable  mass  and  weight  characteristics. 

This  plasma  source  will  be  applied  for  an  ignition  of  hydrocarbon  fuels  in  the  supersonic 

direct-flow  channel  developed  designed  and  manufactured  in  MSU. 
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General  conclusions 

Applicability  of  plasma  generators  with  divergent  plasma  nozzles  for  ignition  of  the 
supersonic  incident  flammable  mixture  flow  has  been  demonstrated. 

Two  embodiments  of  a  stationary  plasma  generator  with  a  vortical  flow  stabilizing  and 
diverging  plasma  ehannel  have  been  realized.  Our  theory  and  experiments  show  that  they  insure 
high  mass  flow  rate,  high  efficieney  of  working  gas  heating  and  smallness  of  heat  losses  into 
water-cooled  anode  surfaee.  The  first  plasma  generator  is  a  large  deviee  for  application  to 
complex  chemically  aggressive  media.  The  second  plasma  generator  is  a  eompaet  apparatus.  It 
can  be  applied  for  ereation  of  a  quasi-plane  jet  in  a  real  seramjet,  using  a  row  of  sueh  plasma 
generators.  Analysis  of  radiation  speetra  of  the  first  generator  plasma  jet  showed  that  the  gas 
temperature  investigated  range  of  inputted  power  in  the  plasma  relaxation  region  (behind  the 
nozzle  entrance)  lies  in  the  range  Tg  =  (3500  -  6600°)  K  sufficient  for  ignition  of  hydrocarbon 
fuels. 

Lagrange  method  on  a  completely  eonservative  implieit  differenee  seheme  with  an  adaptive 
triangular  unstructured  grid  was  used  in  eomputations.  The  total  number  of  knots  was  15000- 
150000.  The  mesh  was  concentrated  in  the  region  of  the  combustion  by  the  adaptive  method. 

This  method  is  applied  for  elarification  of  initial  and  boundary  eonditions  at  investigation  of 
the  plasma  jet-flammable  mix  eross  flow  interaction  (at  angles  135,  150°). 

Computations  of  plasma  jet  injection  to  the  incident  flow  of  propane-air  flammable  mix 
have  been  made.  Plasma  forming  gas  (nitrogen)  mass  flow  rate  was  2  g/s,  plasma  generator 
power  was  15  kW  and  30  kW  at  the  temperature  of  the  ineident  propane-air  stoichiometrie  Our 
ealeulations  show  deep  penetration  of  plasma  jet  into  the  ineident  freestream  and  crossflows  of 
air  and  flammable  mixture.  The  model  of  ignition  and  combustion  of  the  stoichiometrie  propane- 
air  mix  with  a  help  of  the  plasma  jet  from  the  plasma  generator  to  the  ineident  flow  of  the 
flammable  mix  has  been  developed.  It  uses  one  step  effective  chemical  reaction  between 
propane  and  oxygen  moleeules. 

Results  of  computations  demonstrated  the  ignition  possibility  of  propane-air  mix  with  a 
help  of  plasma  generators  with  the  power  1 5  kW  at  the  ineident  propane-air  stoiehiometrie  mix 
temperature  1000  K.  Realizatuion  of  eombustion  wave  towards  the  supersonie  flow  both  in 
eylindrical  and  flat  cases  was  observed  in  computed  variants  at  typical  times  of  ten 
mieroseconds.  The  mixture  is  eompressed  on  the  combustion  front  and  it  in  turn  accelerates 
ignition. 

As  showed  our  previous  computations  concentrations  of  atoms  in  plasma  are  very  high. 
This  can  additionally  increase  the  fuel  activation  in  the  incident  flow. 
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Present  investigations  show  the  principle  possibility  of  supersonic  channel  creation  with 
geometry  of  chamfered  wall  without  a  step,  as  it  is  shown  in  Fig.  Cl.  The  conesponding 
calculations  are  necessary  for  choice  of  optimal  parameters  of  plasma  generator  and  srrpersonic 
channel. 

Calculations  show  that  “high  power  and  low  energy”  generators  can  be  applied  for 
combustion  in  supersonic  conditions. 


Fuel  flow 


Fig.C.l. 


Undertaken  investigations  show  that  appearance  of  0.1%  of  atomic  oxygen  in  the  mix  leads 
to  ~  1.5-2  times  greater  redirction  of  propane  (ignition)  diuing  the  same  time  as  in  the  case  of  the 
imdisturbed  mix.  Appearance  of  1%  of  atomic  oxygen  in  the  mix  leads  to  -5-6  times  greater 
reduction  of  propane  (ignition  becomes  shorter  by  ~5-6)  during  the  same  time  as  in  the  case  of 
the  imdisturbed  mix. 

Moderate  power  self-maintained  discharge  (<10.0  kW)  can  be  used  for  ignition  and 
ignition  enhancement  of  propane-air  stoichiometric  mix  due  to  generation  of  atomic  oxygen  and 
heating  of  the  mix  by  the  discharge.  High  amormt  of  vibrationally  and  electronically  excited 
states  can  considerably  reduce  necessary  power  value.  Application  of  discharges  combination: 
the  non  -selfinaintained  discharge  and  plasma  jet  (see  C.2)  at  low  non-selfmaintained  discharge 
power  (~20  W)  can  be  perspective  for  the  ignition,  since  vibrationally  and  electronically  excited 
oxygen  states  can  significantly  enlarge  reaction  rate  constants  and  can  be  effective  at  interaction 
of  incident  excited  mix  with  the  plasma  jet. 

More  precise  values  of  necessary  ignition  enhancement  power  value  require  detailed 
solution  of  non-selfinaiutaiued  discharge  in  air-firel  mix  including  rates  of  vibrationally  and 
elechonically  excited  states  of  molecular  oxygen  and  hydrocarbon  molecules  reactions  and  the 
solution  of  gasdynamics  equations  in  the  definite  duct.  The  sohrtion  of  the  first  problem  is 
complicated  because  of  necessary  rate  constants  absence.  The  solution  of  the  second  problem 
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can  be  realized  in  approach  developed  in  this  Report  and  can  be  the  subject  of  fijither 
investigations. 


Fig<C2  Schematics  of  the  channel  with  combination  of  discharges 

1  —  supersonic  M=2  nozzle,  2  —  collectors  of  prelimiiiaiy  fiiel  injection  (first  injector  set),  3  -  insulator  channel,  4  - 
collectors  of  main  hiel  injection  (second  injector  set),  5  -  1st  section  of  combustor,  6  —  2  nd  section  of  combustor,  7  - 
quaitz  window,  8  -  third  fuel  injector  set,  9  -  plasma  jet  generator,  10  —  non-selfiiiaintained  discharge  generator 
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